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ABSTRACT 

We measure the net energy efficiency of supernova and stellar wind feedback in the starburst galaxy M82 and 
the degree of mass-loading of the hot gas piston driving its superwind by comparing a large suite of 1 and 2- 
dimensional hydrodynamical models to a set of observational constraints derived from hard X-ray observations 
of the starburst region (the fluxes of the Hea and Lya-like lines of S, Ar, Ca and Fe, along with the total 
diffuse £■ = 2 - 8 keV X-ray luminosity). These are the first direct measurements of the feedback efficiency 
and mass-loading of supernova heated and enriched plasma in a starburst galaxy. We consider a broad range of 
plausible parameters for the M82 starburst, varying the age and mode of star formation, the starburst region size 
and geometry, and supernova metal yields. Over all these varied input parameters all the models that satisfy 
the existing observational constraints have medium to high thermalization efficiencies (30% < e < 100%) 
and the volume-filling wind fluid that flows out of the starburst region is only mildly centrally mass loaded 
(1.0 < /3 < 2.8). These results imply a temperature of the plasma within the starburst region in the range 30 - 
80 million Kelvin, a mass flow rate of the wind fluid out of the starburst region of Mtot ~ 1 .4 - 3.6Mq yr"^ and 
a terminal velocity of the wind in the range Voo = 1410 - 2240 km s"' . This velocity is considerably larger than 
the escape velocity from M82 (vesc ^ 460 km s"' ) and the velocity of the Ha emitting clumps and filaments 
within M82's wind (v'hq ^ 600 km s"' ). Drawing on these results we provide a prescription for implementing 
starburst-driven superwinds in cosmological models of galaxy formation and evolution that more accurately 
represents the energetics of the hot metal-enriched phases than the existing recipes do. 

Subject headings: ISM: bubbles — galaxies: individual : NGC 3034 (M82) — galaxies: halos — (galaxies:) 
intergalactic medium — galaxies: starburst — X-rays: galaxies 



1. INTRODUCTION 

The presence of metals in the intergalactic medium (IGM) 
at all red-shifts ( Songaila 1997). the present-day galaxy mass - 
metallicity relationship ' dGarnettlllOOl iTremonti et alj|2004 . 
and ~ 100-kpc-scale h oles in the IGM at red-shift z ~ 3 
(lAdelberger et alJ[2003h have all been interpreted as the con- 
sequence of energetic metal-bearing outflows from galaxies. 

There is no doubt that at least one form of galactic wind 
exists: starburst galaxies are commo nly observed to host 
supernova-driven galactic-scale winds (iHeckman et al.lfT990l : 
iHec kman 2003) in both the local and h igh redsh ift Universe 
lUehnert & Heck man 1996a; Pettini et al. 2001; S haplev et all 
120031) . These superwinds are phenomenologically complex 
and are comprised of mul tiple different gas pha s es mov- 
ing at different velo cities ("Bland-Hawthorn' ' 1 995'; ' DahlemI 
[T997HHeckman et al . 2000, 2001; Veifleux et al. 2005). How- 
ever, the very hot plasma long hypothesized to be the prime 
mover in superwinds (Chevalier & Clegg 1985), often called 
the wind fluid, has only recently been directly observed in 
the archety pal superwind associated with the starburst g alaxy 
M82 (Grif fiths et aPIIOOOl: Isrickland & Heckmanlllool . 

The most fundamental long-standing uncertainty in our un- 
derstanding of the nature of superwinds and their larger scale 
role is the value of the energy per particle of the metal- 
enriched wind-fluid that drives an individual superwind. Both 
supernovae (SNe) and stellar winds from massive stars are 
responsible for supplying the kinetic energy that drive super- 
winds (Stellar winds only contribute ^ 10% of the total me- 
chanical power of a starburst, but they may supply up to half 
of the gas mass returned to the ISM, see iLeitherer et alJl 19921 ; 

' Department of Physics and Astronomy, The Johns Hopkins University, 
3400 N. Charles St., Baltimore, MD 21218, USA. 



ILeitherer & Heckmanlll995h . What fraction of the --10^' erg 
of kinetic energy released per star that ultimately goes super- 
nova is actually available to drive bulk motions of the ISM? 
As it is understood that SNe dominate the energy return this is 
often termed the SN thermalization efficiency or the SN feed- 
back efficiency, even though stellar winds do play some small 
part energetically. Is the metal-enriched mix of SN and stellar 
wind ejecta diluted and cooled by mixing with the ambient 
ISM, i.e. is the wind-fluid mass-loaded? If this ejecta is to 
escape the galaxy and enter the IGM then its specific energy 
must exceed its initial gravitational potential energy, and fur- 
thermore the total energy in the wind-fluid must be sufficient 
to perform the mechanical work necessary to move the inter- 
stellar and halo gas that lies between the starburst region and 
the IGM (Mac Low & McCray 1988; Strickland et al 2004b). 

For moderate to high thermalization efficiency (efficiencies 
> 10%), the wind-fluid is expected to have an initial temper- 
ature within the starburst re gion i n the range 10 - 100 million 
Kelvin (IChevaher & C legg 1985), even if it has been lightly 
mass-loaded with cold ambient gas. Such temperatures are 
far hotter than the neutral, warm ionized or soft X-ray emit- 
ting plasma s that are the most easily observed phases within 
superwinds dDahlemll 1 997t I VeiUeux et"ai]|2005h . 

However it has only recently become technically possible 
to detect the diffuse hard X-ray emission that would be as- 
sociated with a plasma in the temperature range 10^ < logT 
(K) < 10*^. The starburst regions of even the nearest starburst 
galaxies only subtend ~ 15 - 45" on the sky and are heavily 
populated with X-ray-luminous compact objects. It was only 
with the advent of ^ 1 "-spatial resolution X-ray astronomy in 
1999 following the launc h of the Chandra X-r ay Observatory 
that made it possible for iGriffiths et alJ (120001) to resolve out 
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the point source emission in the archetypal nearby starburst 
galaxy M82 and reveal residual diffuse hard X-ray emission, 
although of uncertain origin. 

It is important to distinguish between the material that ac- 
tually drives a superwind, and the other more easily observed 
phases of swept-up and entrained gas that are part of the wind. 
We define the wind-fluid as the material whose high thermal 
and/or ram pressure drives the superwind. The wind fluid 
certainly comprises a significant fraction of the merged SN 
ejecta and massive star stellar wind material, probably ad- 
mixed (mass-loaded-) with some ambient gas from with star- 
burst region. 

Although other gas phases may dominate the total gas 
mass budget of a superwind, they are less important with 
respect to the issue of energy and metal transfer into the 
IGM. Furthermore the ejection of these phases is less than 
certain, as the measured warm neutral and warm ionized 
medium (WNM, WIM) gas velocities are typically compa- 
rable to the galactic escape veloc ities dHeckman et alj|2000l : 
iMartinI l2005t iRupke et al.] 12005*). Unfortunately even the 
most recent attempts to include superwinds in cosmologi- 
cal models of galaxy formation use only simplified single- 
phase winds, often with energetics and masses based on 
the observed WNM/WIM properties of superwind s (see e.g. 
ISpringel & Hernquisll l2003al: iBertone et all 12005 ). These 
models do not properly address the physics of hot-phase metal 
transport in addition to having overly low values of the energy 
per particle, which casts the reliability of their conclusions 
into doubt. 

Attempts have been made to use measurements of the tem- 
perature of the extended soft X-ray emitting plasmas (T ^ 2 
- 7 X 10^ K) in starburst galaxies to constr ain whether th e hot 
gas will escape the host galaxy (see e.g. iMartinI [19991) . but 
such an approach neglects both the significant kinetic energy 
expected to be associated with this material and that the soft 
X-ray emitting plasma is an indirect probe of the properties of 
the true wind-fluid (IStrickland & Stevensll2000l) . 

The detection of apparently thermal diffuse hard X-ray 
emission of temperature ^ 4 x 10^ K in the starburst region 
of M82 in the 33 ks-long Chandra ACIS-I observations of 
iGriffiths eF al. (2000) suggested it was possible to directly ob- 
serve and constrain the properties of the wind fluid in a star- 
burst superwind. 

We recently obtained and analyzed a new 18 ks Chandra 
ACIS-S observation of M82, and also reprocessed archival 

- We will use the term mass loading to denote the addition and efficient 
mixing of ambient ISM into th e full volume of the wind-fluid, which is con- 
sistent with the terminology in lSuchkov et aljil996 1 Technically mass load- 
ing need not affect the entire flow and can be localized, see e.g. Dyson et al. 
Il993h . Furthermore the form of mass-loading considered in this work oc- 
curs only within the starburst region itself ("centralized mass loading") al- 
though ultimately it affects the entire flow. We find it conceptually beneficial 
to differentiate between mass loading and entrainment in superwinds. Mass 
loading affects the density of one gaseous phase by mixing in material from 
another phase, almost always a cooler and denser phase. The origin and fate 
of the secondary phase is unimportant given the focus on the properties of 
the primary phase. Entrainment is the incorporation into the superwind of 
subsidiary gas phases that leads to creating or sustaining a multiphase super- 
wind. The entrained material is usually ambient disk or halo gas over-run or 
swept-up by the wind fluid, and retains much of its original character, albeit 
with some modification. Although entrainment processes may lead to mass 
loading of the wind fluid the two processes are not synonymous. In princi- 
ple entrainment could occur without leading to significant mass loading of 
the wind fluid. One can not realistically hope to represent or approximate 
the totality of a multiphase superwind as a single phase mass loaded flow. 
Only specific components of the superwind are amenable to such treatment, 
of which the wind fluid is the clearest example. 



Chandra ACIS-I and XMM-Newton observations^ (total ex- 
posure of 48 ks with ACIS-I, 75 ks with XMM-Newton PN 
and 102 ks with each MOS instrument). We confirmed the 
presence of both diffuse continuum and diffuse E ^ 6.7 keV 
Fe Hea line emission within the central starburst region 
(IStrickland & Heckmanll2007i. hence forth referred to as Pa- 
per I). Having assumed that the line-emitting plasma is metal- 
enriched to Zpe ~ 5ZFe,0, we found that the observed iron line 
luminosity of Lx.6.7 (1 . 1 - 1 .7) x 10""^ erg s"' was consistent 
with the properties of the wind fluid expected given M82's 
star formation rate, provided that the SN thermalization effi- 
ciency was high. However the diffuse hard X-ray continuum 
luminosity of Lx.2-8 ~ 4.4 x lO-'^ erg s"' was too luminous to 
be the thermal bremsstrahlung associated with the wind fluid 
given the parameters assumed when interpreting the iron line 
emission. 

In this paper we consider in greater detail what constraints 
the observed hard X-ray properties of M82 can place on the 
properties on the superwind and the wind-fluid, in particu- 
lar the efficiency of SN heating (e) and the degree of central 
mass-loading (/3). Ultimately these determine the energy per 
particle of the wind fluid. We restrict analysis to the star- 
burst region within r ~ 500 pc of the center of M82. We cal- 
culate the X-ray emission using both multi-dimensional hy- 
drodynamical simulations and the 1 -dimensional analytical 
iChevalier & Cl egg (1985) model (hence forth referred to as 
the CC model). By using the simple but relatively accurate 
CC model we can explore the predicted X-ray properties of 
the wind fluid in models covering a wider range of parame- 
ter space than if we used multi-dimensional hydrodynamical 
simulations alone. 

A variety of recent papers have presented theoretical cal- 
culations of the broad-band X-ray emission from individual 
star-forming clusters, using either the adiabatic CC model, 
or the non -adiabatic solution of Silich and collaborators 
dCanto et al 2000; Stevens & Hartwell 2003; Oskinova 20051 
ISilich etal|2005i; iJi et al. 2006), with the aim of seeing if the 
particular model used can explain the observed data. Our mo- 
tivation differs from these papers, in that our aim is to use 
theoretical models in combination with observational data on 
M82 to constrain fundamental properties of the superwind. In 
addition our approach differs from these previous papers in its 
use of X-ray line luminosities rather than purely broad-band 
X-ray luminosities. In reality the broad-band emission might 
be dominated by non-thermal continuum sources. Further- 
more line emission is more sensitive than the thermal contin- 
uum luminosity to interesting parameters such as the plasma 
temperature. 

2. OBSERVATIONAL CONSTRAINTS ON HIGH 
TEMPERATURE GAS IN M82 

The observational constraints derived in Paper I on the dif- 
fuse hard X-ray emission within r = 500 pc of the nucleus of 
M82 from the Chandra and XMM-Newton observations are 
summarized in Table [T] These comprise a robust 6.69 keV 
iron line detection (> 4(t in significance), an upper limit on 
the 6.96 keV iron line, and from the Chandra observations the 
luminosity of the diffuse continuum. For the Chandra-hased 
observations the values refer to the diffuse emission alone, af- 
ter the exclusion of detected point-like X-ray sources and cor- 
rections for undetected point sources based on the observed 

3 Chandra Obsid numbers 2933, 361 and 1302. The XMM-Newton data 
Obslds are 0112290201 and 0206080101. 
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Fig. 1 . — A three color composite image of M82 sliowing tlie region within 
a 500 pc radius of the nucleus (white circle) in comparison to the galaxy and 
superwind. Soft X-ray emission in the 0.3-2.8 keV energy band is shown in 
red, optical R-band emission (starlight) in green, and diffuse hard X-ray emis- 
sion in the 3-7 keV energy band in blue. The X-ray images are adaptively 
smoothed Chandra ACIS-S images with point source emission removed and 
interpolated over. The image is 4f8 {r^ 5 kpc) on a side. 

point source luminosity function. These luminosities also in- 
clude 32% and 20% upward corrections to the luminosities 
derived from the ACIS-S and ACIS-I data to account for re- 
gions excluded from the spectra due to remove the emission 
from detected point sources. For the XMM-Newton observa- 
tion the line luminosities include all emission within r = 500 
pc, but the contribution to the line luminosities from point-like 
sources is expected to be small in this particular observation. 

Note that the region of bright diffuse hard X-ray emission 
studied in this paper is associated with the nuclear starburst re- 
gion. As can be seen from Fig.[T]this region is much smaller 
than the full 12 -kpc-scale superwin d that is most often stud- 
ied in optical (Lehn ert et al.l 11999!) or soft X-ray emission 
(ISteven s et al. 2003). 

As discussed in Paper I the nature of the diffuse broad-band 
continuum is not clear. If the continuum is thermal then the 
relative strength of the iron line emission to the observed con- 
tinuum is much lower than expected for a plasma of Solar 
elemental abundances, let alone the iron-enriched plasma ex- 
pected from recent core-collapse SNe. We demonstrated that 
a wind mass-loaded with additional ambient gas could pro- 
duce the total broad-band X-ray luminosity that we observe, 
but the required degree of mass loading could not sufficiently 
reduce the relative strength of the iron line to the continuum. 
Thus the best fit bremsstrahlung models for the continuum do 
not necessarily provide meaningful estimates of the temper- 
ature of the iron-emitting plasma, and should be interpreted 
with caution. Power law spectral models provide a marginally 
better fit to the continuum than do thermal bremsstrahlung 
models. Thus a non-thermal process such as Inverse Comp- 
ton X-ray emission might be responsible for the majority of 
the diffuse hard X-ray continuum. The X-ray luminosity as- 
sociated with Inverse Compton emission depends sensitively 
on the unknown magnetic field strength within the starburst 
region. In Paper I we estimated it might account for ^25% 
of the observed diffuse X-ray luminosity in the £ = 2 - 8 keV 
energy band, although other authors have in the past given 



higher lu minosity estirnates for Inverse Compton em ission 
from M82 (ISchaaf et al.lfT98llMoran & Lehnertlll997h . 

The upper limits on the iron 6.96 keV / 6.69 keV line flux 
ratio"^ of < 0.34 place upper limits on the temperature of 
T < 7.8 X lO^K if the plasma is in collisional equilibrium 
(see Table [2j. The Coulomb relaxation and ionization equi- 
librium timescales presented in Paper I are much shorter than 
the timescale for plasma to flow out of the starburst region, 
so the assumption of collisional ionization equilibrium is rea- 
sonable. 

Other lines in the 2<E (keV) < 10 hard X-ray band may in 
principle be used to constrain the plasma temperature in a star- 
burst region. In particular the helium-like and hydrogen-like 
ions of sulphur (at £ ^ 2.45 keV and E ^ 2.62 keV respec- 
tively), argon (£ 3.13 keV and £ 3.33 keV) and calcium 
(E - 3.90 keV and £ - 4.11 keV) can be detected in CCD 
specti-a from the Chandra ACIS and XMM-Newton EPIC in- 
struments. These lines are most intense in plasmas with tem- 
peratures in the range logT ^ 7.0 - 7.5, somewhat lower than 
the optimal temperature for the iron lines, but this might be 
advantageous in diagnosing starbursts where mass loading is 
significant and/or the efficiency of SN thermalization is low. 

In Paper I we specifically excluded data at the expected en- 
ergies of the S, Ar and Ca lines from the fitting of the M82 
nuclear region spectra. Repeating the analysis described Pa- 
per I but now fitting all the data in the energy range E = 2.2- 
9.0 keV in both the XMM-Newton and Chandra ACIS-S spec- 
tra yields the line detection significances and line luminosities 
given in Table |3] We fit for the intensity of the helium-like and 
hydrogen-like S, Ar and Ca lines, having fixed the line ener- 
gies at the expected values. These lines are faint, but in some 
cases likely to be real. As an example we display the best-fit 
to the XMM-Newton EPIC pn H- MOS spectra in Fig. |2] De- 
tection significance was assessed using Monte Carlo methods 
with several thousand simulated spectra. Where the line de- 
tection significance was < 3cr we give the 99% upper limit on 
the line luminosity (and not the best-fit value). The Chandra 
ACIS-I spectra gave similar results to those from ACIS-S, but 
as the line detections are less significant than in the XMM- 
Newton data there is little reason to present the values. 

In the XMM-Newton EPIC pn + MOS spectra of M82's nu- 
cleus the equivalent of the Hea line from highly ionized S, Ar, 
Ca and Fe are all detected with significances > 3<t, while in 
the shorter (and hence lower S/N) Chandra ACIS-S spectrum 
of the nuclear diffuse emission the S Hea and Fe Hea are de- 
tected at > 3(7 while the Ar Hea and Ca Hea lines are only 
detected at < 2a. In both instruments Si Ly/3 is also detected 
at high significance, but we shall ignore this lower ionization 
state line in the analysis to follow. The only Lya-like line 
detected at > 3ct is S Lya in the XMM-Newton EPIC specti-a. 

Using the robust line detections (i.e. those with significance 
> 3(7) we can estimate the plasma temperature from the ratio 
of S Lya to S Hea, and from between different elements from 
the ratio of the various Hea lines (assuming their relative el- 
emental abundances are Solar). These temperature estimates 
are presented in Table |2] The S ions tightly constrain the ion 
temperature to logT w 7.1, while the ratio of the S, Ar or Ca 
Hea line to the Fe Hea line yields logT ~ 7.3. The ratio 
of S Hea to either Ar Hea or Ca Hea yields a temperature 
estimate of log T ^1 .2. 

This pattern of the inferred temperature increasing as higher 

■* Note that the line ratios quoted in Paper I were the ratio of the photon 
fluxes, whereas in this paper we quote the energy flux ratio. 
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TABLE 1 

Diffuse continuum and iron line luminosities from Paper 1 . 



Instrument 
(1) 


ObsID 

(2) 


logL2_8keV 
(3) 


logLpcKa 

(4) 


(5) 


log^FeLya 

(6) 


XMM-Newton EPIC 
Chandra ACIS-S 
Chandra ACIS-I 


0206080101 
2933 
361+1302 


-,o £5+0.02 
39 64+0.02 


37.31+0-23 
37 g4+8:S 
J '•^'+-0.48 
37 03+0.20 
^'•"-0.26 


38.23+[j-^6 
38.13+0;^^ 


< 37.66 

< 38.15 

< 38.01 



Note. — This table summarizes the results presented in Strickland & Heckman 1 2007) on the emission 
from within r = 500 pc of the nucleus of M82. Column 2: The ObsID is the identification number allocated 
to each observation by the Chandra and XMM-Newton Science Centers. Column 3: Logarithm of the total 
diffuse X-ray luminosity in the £ = 2 - 8 keV energy band, including line emission. Column 4: Logarithm 
of the E ~ 6.4 keV Fe Kq line luminosity. Column 5: Logarithm of the E ~ 6.7 keV Fe Hea line 
luminosity. Column 6: Logarithm of the E ~ 6.96 keV Fe Lya line luminosity. The quoted uncertainties 
are 68.3% confidence for one interesting parameter, while the upper limits are 99.0% confidence. 



TABLE 2 

Temperature estimates for the diffuse hard X-ray emission. 



Basis for temperature estimate 


Detector 


logr(K) 


(1) 


(2) 


(3) 



Fe Lya / Fe Hea 


XMM-Newton 


<7.89 


Fe Lya / Fe Hea 


ACIS-S 


< 8.05 


Fe Lya / Fe Hea 


ACIS-I 


< 8.04 


Ca Hea / Fe Hea 


XMM-Newton 


7-11+0.10 
'•-'^-0.06 


Ca Hea / Fe Hea 


ACIS-S 


>7.21 


Ar Hea / Fe Hea 


XMM-Newton 


7.30 ±0.03 


Ar Hea / Fe Hea 


ACIS-S 


>7.29 


S Hea / Fe Hea 


XMM-Newton 


7.26 ±0.01 


S Hea / Fe Hea 


ACIS-S 


7 9^+0.06 
'■■^-^-0.03 


(Ca Lya + Ca Hea) / Fe Hea 


XMM-Newton 


> 7.22 


(Ca Lya + Ca Hea) / Fe Hea 


ACIS-S 


> 7.12 


(Ar Lya + Ar Hea) / Fe Hea 


XMM-Newton 


>7.26 


(Ar Lya + Ar Hea) / Fe Hea 


ACIS-S 


>7.20 


(S Lya + S Hea) / Fe Hea 


XMM-Newton 


7 30+° °' 

'•-^"-0.02 


(S Lya + S Hea) / Fe Hea 


ACIS-S 


>7.26 


S Hea / Ca Hea 


XMM-Newton 


7 o-i+O.O? 
' -^^-0.05 


S Hea / Ca Hea 


ACIS-S 


<7.27 


S Hea / Ar Hea 


XMM-Newton 


7. 19 ±0.05 


S Hea / Ar Hea 


ACIS-S 


<7.11 


Ar Hea / Ca Hea 


XMM-Newton 


7.28!0;13 


Ca Lya / Ca Hea 


XMM-Newton 


< 7.63 


Ca Lya / Ca Hea 


ACIS-S 


<7.69 


Ar Lya / Ar Hea 


XMM-Newton 


<7.42 


Ar Lya / Ar Hea 


ACIS-S 


<7.59 


S Lya / S Hea 


XMM-Newton 


7 10+0'°2 


S Lya / S Hea 


ACIS-S 


<7.12 


Bremsstrahlung fit to diffuse continuum 


ACIS-S 


7 64+<'-22 


Bremsstrahlung fit to diffuse continuum 


ACIS-I 


7 50+°'°' 

-0.07 



Note. — Temperature estimates are given only for line ratios where one or more 
of the lines involved is detected at > 3cr. The quoted uncertainties are 68.3% con- 
fidence for one interesting parameter, while both upper and lower limits are 99.0% 
confidence. 



ionization state ions are used (which are of course more sensi- 
tive to higher temperatures) indicates that either the plasma is 
muhi-phase, and/or the the plasma is not in collisional ion- 
ization equilibrium (CIE). Of these explanations we favor 
the former, as multi-dimensional numerical hydrodynamical 
simulations predict a broad range of gas temperatures below 
logT ^ 8 (see e.g. [Strickland & Steven s 2000). Furthermore 
we discussed in Paper I reasons to expect the Fe-emitting 
plasma to be in or close to CIE. 

Unfortunately these results do not strongly constrain the 
temperature of the Fe-emitting plasma. One possibility is that 
the iron line emission arises in the very hottest phases in the 
superwind, i.e. at logT < 7.9, and that most of the S, Ar and 



Ca emission comes from somewhat cooler components asso- 
ciated with the wind/ISM interaction and the larger scale soft 
X-ray emission in the M82 superwind. In this case, given the 
possible soft X-ray contribution, the S, Ar and Ca observed 
luminosities must be treated as upper limits to the S, Ar and 
Ca line emission coming from the volume filling wind-fluid. 
Alternatively the wind fluid is relatively cool, with a tempera- 
ture corresponding to that measured from the Ar and Ca lines, 
i.e. \ogT ^ 7.3, and the Fe lines are produced by some unex- 
plained non-thermal mechanism. 

In either case a successful model of M82's nuclear emission 
must reproduce the observed hard X-ray iron line emission 
without over-producing the broad-band continuum and S, Ar 
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Fig. 2. — (Upper panel) The XMM-Newton EPIC pn (black data points with error bars), MOSl (red data points) and MOS2 (green data points) spectra of the 
central 500 pc of M82. The best-fit spectral model (red solid line) is shown, along with the named X-ray line components of the model (dashed and dotted fines 
with a color matching the parent data). The Hea-equivalent lines of S, Ai', Ca and Fe, along with Si Ly/3 and S Lya, are all detected with significances of > 3(T. 
The other lines have lower significance. (Lower panel) The residuals from the best-fit model. 



TABLE 3 

Hard X-ray S, Ar and Ca line detections and luminosities in M82. 



Line 
(1) 



Ion 

(2) 



Energy 
(3) 



Detector 
(4) 



Significance 
(5) 



logLx 
(6) 



Hea 
SLya 
Ar Hea 
Ar Lya 
CaHea 
Ca Lya 
SHea 
SLya 
Ar Hea 
Ar Ly« 
CaHea 
Ca Lya 



S XV 
SXVI 
Ar XVII 
Ar XVIII 
Ca XIX 
CaXX 

sxv 

SXVI 
Ar XVII 
Ar XVIII 
Ca XIX 
CaXX 



2.45 
2.62 
3.13 
3.32 
3.90 
4.11 
2.45 
2.62 
3.13 
3.32 
3.90 
4.11 



XMM-Newton 

XMM-Newton 

XMM-Newton 

XMM-Newton 

XMM-Newton 

XMM-Newton 

ACIS-S 

ACIS-S 

ACIS-S 

ACIS-S 

ACIS-S 

ACIS-S 



> 3.5 

>3.5 

>3.5 

0.8 

3.5 

1.1 

>3.9 

1.7 

1.9 

0.5 

2.4 

0.5 



38.19±0.02 

->-, 43+0.08 

37 53+"-*" 
-"■•^ -0.08 

< 37.26 

9-7+0.12 

^'■^'-0.17 

< 37.14 
38.32+|;-5;« 

< 37.64 

< 37.56 

< 37.69 

< 37.52 

< 37.54 



Note. — The luminosities are given as logjg of the line luminosity in units of ergs" , 
assuming a distance to M82 of Z) = 3.6 Mpc and zero intrinsic absorption. Line energies are 
in units of keV. Line significances are the cr values equivalent to the enclosed probability \ - p 
for a two-tailed Gaussian, where p is the probability of the line being spurious. In the case of 
highly significant lines we can only estimate a lower limit to the significance of the line. The 
quoted uncertainties on the luminosities are 68.3% confidence, while upper limits on logLx 
are provided at 99.0% confidence if the line significance is < 3o'. 



and Ca line emission. 

3. THEORETICAL AND NUMERICAL 
METHODOLOGY 

Hard X-ray emission in multi-dimensional hydrodynami- 
cal models of superwinds is dominated b y the thermal emis- 
sion fr om starburst region itself (Tomisaka & Bre gmanll 19931 : 
ISuchkov et al. 1994: Strickland & Stevens 2000) . 

If the starburst region is spherical then the plasma proper- 



ties of this region obtained from hydrodynamical simulations 
agree very well with the simple one-dimensional radially- 
symmetric analytical CC model, in which the mechanical 
energy and returned gaseous ejecta from core-collapse su- 
pernovae are injected at a uniform rate only within a radius 
and radiative energy losses are assumed to be negligible. 
Gravity and any surrounding ambient gaseous medium are ig- 
nored. 
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As the feedback parameters we wish to constrain are related 
to the observable X-ray emission via the plasma properties in 
the starburst region it is worthwhile for us to expl icitly repro- 
duce t he analytical solution obtained by Cheval ier & Clegg 
([1985). We consider the conditions under which the adiabatic 
assumption inherent in the CC model may be invalid in Ap- 
pendix [A) but we will show that for most reasonable param- 
eters associated with the M82 starburst (described below) the 
radiative energy losses are not significant. 

In real galaxies starburst regions are better described as 
disks or ring of star-formation rather than spherical regions. 
To explore the predicted X-ray emission for our default 
model of the M82's disk-like starburst we begin by using 2- 
dimensional hydrodynamical simulations in cylindrical sym- 
metry to predict the plasma properties within the a radius of 
500 pc of the center of the galaxy. 

However, with appropriate scaling (see Appendix |B]i, it is 
possible to use the 1 -dimensional CC model with a spheri- 
cal starburst to predict the plasma properties of the wind fluid 
within a non-spherical starburst region to high accuracy, and 
to predict the hard X-ray emission from even larger regions 
with an accuracy no worse than the uncertainties in the ob- 
served X-ray luminosities. 

These scaled CC models are considerably faster to calcu- 
late than a full hydrodynamical model, so we will use them to 
explore how the constraints on the thermalization efficiency 
and mass loading change if we adopt a different range of pa- 
rameters for the M82 starburst. 

This approach, in particular the use of the scaled CC model, 
is applicable only to the study of the wind fluid in the vicinity 
of the starburst region. It is not intended for use as a model 
of an entire superwind as it neglects the interaction between 
the wind fluid and the multiphase ambient ISM. The simplifi- 
cation of neglecting the ambient ISM works because the only 
significant source of thermal hard X-ray emission in a super- 
wind is the wind fluid itself, in particular the wind fluid within 
the starburst region. The interactions between the wind fluid 
and the multiphase ISM are important in the large scale wind 
and for its emission at soft X-ray and longer w avelengths (e.g. 
[Strickland & Stevensll2000t ICooper et al.ll200 8). but with the 
exception of altering the net efficiency of SNh-S W heating and 
any mass-loading of the wind fluid these interactions are not 
expected to alter the properties of the wind fluid within the 
starburst region to any significant degree given its energetic 
dominance. 

3.1. 2-Dimensional Hydrodynamical Modeling 

The hydrodynamical simulations were performed using 
th e multi-dimensional numerical hydrodynamics code VH- 
1 (iBlondin et al.lll99 ^. Simulations were performed in 2- 
dimensional cylindrical symmetry with a cell size of 2 pc 
over 500 by 500 cell computational grid. Reflecting bound- 
ary conditions were imposed along the r and z axes, with in- 
flow/outflow boundary conditions are the other edges of the 
computational grid. 

Energy and mass are inj ected within the starburst re gion 
in the manner described in [Strickland & StevensI (l2000l) . In 
general the starburst region was modeled as disk-like region of 
the size described in § 13.31 No gravitational forces or radiative 
cooling are considered for the purposes of consistency with 
the analytical CC model described below. 

The computational grid was initially populated with a ten- 
uous ISM of density p = 2 x 10"^^ gem and temperature 
r = 8 X 10^ K. The simulations were run until the wind had 



blown the initial ISM off the grid and reached a steady-state 
solution, which typically occurred f ^ 1 - 2 Myr after the on- 
set of the simulation. Gas properties and X-ray luminosities 
are measured at a time of 3 Myr within a spherical region of 
radius of 500 pc to closely match the observed region dis- 
cussed in § ID 

3.2. The 1 -Dimensional Analytical CC Solution 

For completeness we present the analytical flow solution 
derived by CC, and how we use it to derive the fluid proper- 
ties as a function of radial position. More complicated, but 
also more g eneralized, method s of deriving the solutio n are 
presented in ICanto et all (l2000h and ISihch et al.1 (|2004|) . The 
basic model is completely specified by only three parameters: 
the energy and mass injection rates (ZstotjA^tot); and the radius 
of the starburst region {Ri,) within which the energy and mass 
are injected. 

CC present the stationary solution to the hydrodynamical 
equations as a function of the Mach number M = v/cj, where 
V is the local speed of the flow and the sound speed is Cs = 
y/i'yP/p). The pressure P, density p and adiabatic index 7 
have their normal meanings. At radii r < mass and energy 
are injected at a uniform constant rate per unit volume qnj = 
Mtot/Vi, and qE = -Etot/V*, where the starburst volume = 
47r/?^/3, and the solution for the Mach number is 



37+1/M2 



(37+1) 
(57+1) 



7-1+2/M2 



1+37 J V 1+7 

This approach to the energy and mass injection im- 
plicitly assumes that the exact manner in which stars re- 
turn mechanical energy to the ISM can be ignored due to 
rapid local thermalization and mixing, i.e. individual stel- 
lar winds and SN blast waves collide and merge on spa- 
tial scales small compared to the size of the starburst re- 
gion. If radiative energy losses are to be small then these 
collisions must occur with a time shorter than the cool- 
ing time scale for stellar wind bubbles or supernova rem- 
nants (SNRs). Theore tical inv estigations of some of these is - 
sues may be found in Larson ( 1974), Heckman et al. (1990^, 
ICantoetalJ (l20()0h. .M eHoU & de GouveiaJDaLPino (200i) 
and lMelioli et al.1 ( l2005h . 

At larger radii r > then qM = qE = 0, and the solution is 



7-1+2/M- 
1+7 



2\ 2(7-1) 



(2) 



For any given fractional radius r/R^ the Mach number is 
found using a numerical root-finding technique (in our case 
the Newton-Raphson method) on the appropriate equation. 

To obtain the density at any radius r < we integrate 
Equation 1 of CC to obtain 

qMr 

P=-^- (3) 

3 V 

Within the starburst region the gas density, pressure and 
temperature vary only slowly with radius, with maxima at 
r = 0. The central density is effectively determined by the rate 
at which the injected material can flow out of the starburst 
region, cx Mtot/(47r/?^ Vqo) where v^c is the terminal veloc- 
ity of the wind at large radius, and the central temperature is 
set by the energy injected per particle. More specifically, the 
central temperature and density are Tc = QA{pmi{Etot)/{kMtot) 
and Pc = 0.296 Mttf£;'/^7?;2 for 7 = 5/3. 
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At radii r > R^, the mass flow rate through the spherical 
surface of radius r is equal to the total mass injection rate 
within the starburst region, so the density is 



(4) 



To derive the local velocity of the flow v we must first deter- 
mine the speed of sound. If we make the assumption that the 
flow is adiabatic then we can obtain the speed of sound at any 
point (r < Rj, or r > R^,), as we know the Mach number and 
the initial energy per particle is conserved. Specifically we 
solve Equation 3 of CC and substitute for p using Equation[3] 
to obtain 

'^.^)=4^^. (5) 

The speed of the flow at any point is then simply v = Mcj. 
For r < R^, the flow is subsonic, with m cx r. The sonic point 
is reached at r = R^,, beyond which the flow velocity rapidly 

•1/2 '-1/2 

tends toward the terminal velocity Vqo = y2£'tot ^tot ■ 

With the density p and the speed of sound Cj determined 
at every point within the flow the thermal pressure P is also 
now known. The kinetic temperature T = fim-aP/ pk, where 
finiH ~ 1.02 X lO'^'^gcm"^ is the mean mass per particle for 
a fully ionized plasma of Solar composition. 

Numerically we derive the gas density, temperature and ve- 
locity as a function of radius, iterating outward in thin radial 
shells. The numerical implementation we use is based on 
code originally created by Ian Stevens jStevens & Hartwell 
l2003h . This semi-analytical code reproduces the density and 
pressure to within 0.1% of the values predicted by the full 
2-D hydrodynamical simulation using the VH-1 code for a 
spherical starburst region. The scaling used to approximate a 
non-spherical starburst region using the spherical CC model 
is described in AppendixiSl 

Figure [3] presents the radial variation in several fluid vari- 
ables and associated derived quantities for one particular set 
of input parameters, as predicted using our 1 -dimensional CC 
model code. 

3.3. Parameters of the M82 nuclear starburst 

3.3.1. The starburst age and magnitude 

The flow times for the region of interest (r < 500 pc) are 
short (a few Myr) compared to the lifetime of the massive stars 
powering the flow (~ 30 Myr). The properties of the starburst 
region will very rapidly approach those of the steady-state CC 
solution after any change in the energy and mass injection 
rate. 

Nevertheless, in order to apply the CC model to M82 we 
must determine the history (i.e. age and star formation rate as 
a function of time) of the starburst event in order to derive the 
current mechanical energy and mass injection rates and the 
elemental abundances in the injected gas. Furthermore, de- 
termining the degree to which the M82 starburst event can be 
treated as one the simple limiting cases of either a single in- 
stantaneous burst (SIB) or as continuous star formation (CSF) 
is imp ortant, in particular for comparison to Starburst99 mod- 
els (.Leitherer etal.ll 19991) . 

Recent observational studies have greatly expanded our 
knowledge of M82's global star formation history since its 
close enco unter with M 81 approximately 0.6 — 1 Gyr ago 
(see e.g. de Grijs' 2001). Nevertheless significant gaps still 
remain in our understanding of the history of the nuclear star- 
burst. 



Exterior to the central kiloparsec star formation appears to 
have peaked ^ 600 Myr ago in a burst of comparable intensity 
to the current nuclear starburst event, but star fo rmation effec- 
tively ceased in the disk 20 - 30 Myr ago dde Griis et alj 
'2001; Mayya et al. 2006). 

Within the central kiloparsec (projected separations from 
the nucleus of r < 500 pc) the majority of recent studies 
have concentrated on the ^ 200 compact star clusters found 
there dO'Connell et al.l[T99l iMelo et al.ll2005h . The bright- 
est 20 Star cl usters typi cally ap pear to ha ve ages ^5 — 
10 Myr (Satvapal et al.l ll997; Mc Cradv et al.ir2003) . although 
two massive clusters (M82-F and M82-L) at the edge of this 
region h ave ages o f 60 ±20 Mvr (ISmith & G anagheij|200ll: 
ISmith e t al. 20061). iForster Schreiber et alj J2003a b) present 
a detailed analysis of the NIR emission (1 — 4-5 pm) from a 
variety of regions within r < 250 pc, and find that two bursts 
(of ages ~ 5 and ~ 10 Myr respectively) are required in each 
region to simultaneously explain the strength of both nebular 
emission features and CO absorption from Red Super Giant 
(RSG) stars. This work, and the estimated ages of the bright- 
est star clusters, indicates that the nuclear star formation is not 
well represented as a single instantaneous starburst event. 

Note that the cluster-based studies do not provide a com- 
plete view of the starburst, as they are mainly observed in 
regions of lo wer-than-average extinction^ toward the nucleus 
(|Smith et al . 2006a). Nor does this work strongly constrain 
the star formation rate between 10 — 60 Myr ago. Old clus- 
ters in this age range, of similar mass to the bright clusters, 
would be considerably fainter and thus harder to detect than 
the f ^ 10 Myr-old clusters'', due to the rapid decrease of the 
ionizing and bolometric output at ages t >10 Myr UV spec- 
troscopic studies of nearby low-obscuration starbursts have 
shown that the massive star population within the clusters is 
typic ally younger the the field massive star p opulation (see 
e.g. ITremonti etalJl2001l: IChandar et ai]|200 5). The proba- 
ble cause for this effect is the gradual dissolution of the mas- 
sive star clusters on times scales of "^10 — - 30 Myr (for a 
recent review see lde Griis & Parmentieiil2007L and references 
therein). We are not aware of any similar UV study of the field 
stellar population in M82's nucleus, so there is no direct ev- 
idence for the presence of star formation intermediate in age 
between the f ^ 10 Myr-old clusters and the t ^ 60 Myr-old 
clusters. Nevertheless it remains very plausible that signifi- 
cant star formation continued in between the two periods. 

Another reason to suspect that the nuclear burst has 
been active for longer than ^10 Myr is the dynami- 
cal age of the Ha cap, a large nebular feature associ- 
ated with the superwin d at a height z ^ 12 kpc above the 
mid-p lane of the M82 dPevine & Ball vl 1 19991: iLehnert et af] 
Il999l) . Its true space velocity is unknown, but assuming 
a mean velocity of VHa = 600 km s"' equivalent to that of 
the Ha filaments in the lower w ind ( McK eith et al.l 119951 : 
IShopbell & Bland-Hawthoriill 19981) it would have taken ~ 20 
Myr to reach its current location if its origin was within the 
disk of the galaxy. 

^ The typical measured extinction towards the clusters is Ay ~ 2 — 6 
magnitudes, see Melo et al. (2005). 

* Note that this statement, and the previously mentioned cluster ages, is 
based on single star evolutionary models. Mas-Hesse & Cervifio 1 1999) ar- 
gue that if binary systems are taken into account the strongly ionizing phase 
of a single coeval cluster can last up to 30 Myr after formation, and that 
the cluster can simultaneously produce nebular emission and have a substan- 
tial number of RSGs. If this is indeed the case then the clusters identified as 
having ages of ~ 10 Myr may be as old as ~ 30 Myr 
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R (pc) R (pc) R (pc) 

Fig. 3. — (a) An example of the radial CC solution described in § 13.21 Note the relatively uniform number density, thermal pressure and temperature within the 
assumed starburst region of radius = 300 pc, and rapid drop in those parameters outside the starburst region. The net energy and mass injection rates assumed 
are Etot = 3.1 X 10*^ ergs"' and Miot = 1.4Mq yr"'. (b) Local wind velocity v and sound speed Cs. (c) Logarithms of the ionization-state-related parameter r] 
(Equation ll4t and cumulative volume emission integral EI(< R). The plasma is in ionization equilibrium when logr; > 12. Values are given in c.g.s units unless 
it is stated otherwise. 



We thus choose to treat M82's nuclear starburst as having 
experienced a uniform rate of star formation since the onset of 
the burst fburst years ago, with negligible star formation before 
that time. We will treat fbmst as a model parameter in the range 
10 — 50 Myr. 

In Table |4] we present the bolometric luminosity, and total 
SN plus stellar wind energy and mass injection rates (all per 
unit star formation rate) for continuous star formation models 
calculated using ver sion 5.1 of S tarburst99 population syn- 
thesis code (Leithere r et al.lfT999l) . The star formation rate is 
5=1 Mq yr"' for a S alpeter I MF between the mass limits of 
1 — IOOMq. For a'Kroupa! ( 2001) IMF between the mass 
limits 0.01 — IOOM0 producing the same total mass of mas- 
sive stars (mass range 8 — 100 M0) as in the Salpeter IMF 
model, the star formation rate is a factor 1.76 times larger 
The Geneva stellar evolutionary tracks for an initial metallic- 
ity equal to Solar, and with high stellar wind mass loss rates, 
were used. 

The total 8 — 1000//m IR luminosity of M82 is Lir = 
5.8 X 10'" based on the latest re-calibrated 12 - lOO^m 
IRAS fluxes (ISanders et al.ll2003h . IR emission dominates 
M82's spectral energy distribution, so by equating Lir to the 
starburst bolometric luminosity for each assumed age yields 
star formation rates S between 4 — 6 Mq yr"' (Salpeter IMF). 
These values and the associated SN plus stellar wind energy 
and mass injection rates are ta bulated in Table |4] By way 
of comparison the formulae in iKennicutl ([T998), which as- 
sume CSF with a large age, give a star formation rate of 
S = 3.9 Mq yr"' for M82 based on Lir, when corrected to the 
1 — 100 M0 mass range Salpeter IMF we use. 

These star formation rates are not unreasonable. Over the 
time scale of interest the net mass of stars formed is at maxi- 
mum ~ 50 % of the dynami cal mass of ~ 7 x 10^ Mq within 
r = 250 pc dGotz et all 1 199011) . even if we assume a Kroupa 
IMF extending down to a lower mass limit of 0.01 Mq. 

3.3.2. The total size of the starburst region 

The brightest region of diffuse non-thermal 5 GHz radio 
emission in M82 subtends an angular extent of ~ 43" x 
6" (~ 750pc X 105 pc), the same region occupied by the 
compact radio sources presumed to be you ng SN remnants 
jKronberg et al.l 119851; iMuxlow et all 11994 . This is simi- 
lar to the size of the region occupied by the brightest star 
clusters seen in the N IR, - 58" x 8" (-- 1015pc x 140pc, 
iMcCradv etani200l . These major-axis angular diameters 
closely match the base of the superwind as measured by the 



region of spht optical em ission Unes, which has a diameter of 
~50" lG5tzetal]([T990l) . 

If we approximate the starburst region as a flat cylinder 
of diameter di, and total height /z* then the starburst volume 

= 7r(t/*/2)^/!* and surface area 5* = 7rt/*(/!* + t/*/2) in each 
case is y,,i = 4.6 x 10^ pc^ = 1.1 x 10^ pc^ (5 GHz) and 
y*,2 = 1.1 X lO** pc^ , 5*,2 = 2.1 X 10^ pc2 (NIR). The equiv- 
alent spherical CC model should then have a starburst region 
radius of 7?* i = 300 pc or Ri,2 = 405 pc (AppendixlBb. 

The starburst in M82 is p ossibly concentrated along 
the stellar bar (see e .g. lAchtermann & Lacvl Il995t 
IWestmoquette et ani2007l) . and hence is more geometrically 
complex than a uniform disk. Furthermore star formation may 
not have occurred simultaneously throughout the full starburst 
region, but evidence for inwardly or outwardly propagating 
star formati on within r < 500 pc of the n ucleus remains in- 
conclusive ( iForster Schreiber et al.ll2003a h. Exploring the full 
significance of these effects requires 3-dimensional hydrody- 
namical simulations that are outside the scope of this paper. 

However we can investigate the significance of concen- 
trating the supernova activity within a smaller disk-like re- 
gion than the two considered above. We ad opt the starburst 
disk geometry used in Strickla nd & StevensI (|2000) which has 
d^, = 300 pc and hj, = 60 pc. This produces significantly lower 
starburst volumes and surface areas: = 4.2 X 10'' pc^ and 
S^,3 = 2.0 X lO'' pc^. The equivalent spherical starburst radius 
is 3 = 126 pc. 

For the majority of the models we consider we adopt the 
first starburst region geometry described above: a disk of di- 
ameter dj, = 750 pc and total height = 105 pc (or i = 300 
pc for the scaled CC models). The effect of the size of the 
starburst region is explored using the scaled CC model with 
radii 2 = 405 pc and Ri^ j, = 126 pc. 

3.3.3. The effective energy and mass injection rates 

The most significant contributors to the energy and mass 
return rates in a starburst event are core-collapse SNe and 
the stellar winds from the massive stars. Stellar winds are 
thought to contribute an order of magnitude less mechani- 
cal energy than the SNe (assuming equal thermalization ef- 
ficiencies), but they may contribute up to ^ 50% of the total 
mass returned to the ISM ( Leitherer. Robert. & Drissenll 19921: 
Leitherer & Heckman 1995|). As described above, we used 
the Starburst99 population synthesis code to find the com- 
bined SNe plus stellar wind energy and mass injections rates 
for a continuous burst of star formation as a function of age. 



Supernova feedback efficiency and mass loading in M82 9 



TABLE 4 

Continuous star formation models for M82. 



Age 


Lboi/S 


£'sN+sw/5 


*?SN+sw/5 


S 


£'SN+SW 




(Myr) 


(LqMq-i yr) 


( erg s"' Mq^' yr) 


( ■■■ ) 


(Meyr-') 


( erg s"' ) 


(Mq yr-') 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


10 


9.6 X 10' 


3.1 X 10*' 


0.19 


6.1 


1.9 X 10*2 


1.1 


20 


1.1 X 10'" 


5.1 X 10*' 


0.26 


5.1 


2.6 X 10*2 


1.3 


30 


1.2 X lO'" 


6.5 X 10*' 


0.29 


4.7 


3.1 X 10*2 


1.4 


40 


1.3 X 10'" 


7.1 X 10*' 


0.31 


4.5 


3.2 X 10*2 


1.4 


50 


1.3 X 10'" 


6.9 X 10*' 


0.31 


4.4 


3.0 X 10*2 


1.3 



Note. — Column 1 : Age since the onset of continuous star formation. Column 2: Starburst radiated bolometric 
luminosity per unit star formation rate. Column 3: Mechanical energy injection rate due to SNe and stellar winds per 
unit star formation rate. Coluinn 4: Mass injection rate due to SNe and stellar winds per unit star formation rate. Note 
that the values in Columns 2-4 are specific to a Salpeter IMF between the lower and upper inass limits of 1 — 1 00 Mq . 
Column 5: Star formation rate <S derived from L|R/(Lboi/*S), where Lir is the observed IR luminosity of M82. Column 
6: The mechanical energy injection rate associated with this SFR. Column 7: The mass injection rate associated with 
this SFR. See § |3.3.I| for details on the stellar evolutionary tracks and stellar wind models used. 



Note that the values of /ssn+sw and Msn+sw given in Table |4] 
have not been corrected for any initial mechanical energy 
losses due to radiation before thermalization occurs (the ther- 
malization efficiency e) or for mass loading by the addition of 
cold ambient ISM (with a rate Mcoid)- 

The thermalization efficiency and degree of mass loading 
are expected to be a function of the local gaseous environment 
the massive stars find themselves in, as well as the local su- 
pernova rate per unit volume and the density of massive stars. 
Low ambient gas density and/or high SN rate per unit volume 
should lead to a higher hot gas filling factor, which in turns 
increases the thermalizati on efficiency by redu cing the initial 
radiative energy losses ( Heckman et al.lll990 ). However, at 
very high stellar densities the SN and stellar wind ejecta may 
itself become dense enough to cause drastic energ y losses, e.g. 
within the central regions of su per star clusters dSilich etal] 
I2004t iTenorio-Tagle et al.ll2007h . 

To account for these factors we assume that only a fraction 
( of all the SNe and stellar winds occur in regions where con- 
ditions lead to non-negligible thermalization, and the mean 
thermalization efficiency in these regions is what we define as 
e. Essentially <^ is a modification to the total star formation 
rate of the starburst to account for localized spatial variations 
within the starburst of the SN and stellar wind feedback. We 
will refer to ( as the participation fraction, as it is effectively 
the fraction of the total star formation that is actually involved 
in driving the superwind. For the remaining fraction 1-C of the 
SNe and stellar winds we assume that radiative energy losses 
are so significant that no mechanical energy is supplied to the 
starburst region, and that the ejecta never becomes part of the 
outflow. We treat ^ as a model variable, and refer the reader 
to Appendix|C]for a discussion of plausible values for ( given 
the gas content of the starburst and the star cluster mass func- 
tion. 

The net energy injection rate into the hot volume-filling gas 
in the starburst region is 



E tot - ^CE sn+sw • 
The net mass injection rate is 



M,ot = CM: 



sn+sw 



+ Mcold = /?CA^SN+SW, 



(6) 



(7) 



where the mass loading factor /3 is a parameter commonly 
used to indicate the relative amount of mass loading. If /3 = 1 
then the wind is not mass loaded. In this paper we con- 
sider only mass loading that operates within the starburst re- 



gion, i.e. centr al mass loa ding in the terminology adop ted by 
ISuchkov et all ( Il996h and I Strickland & StevenslfcOOOt) . 

Written as a function of the thermalization efficiency, the 
mass loading factor and participation factor, the central tem- 
perature, pressure and density are 



7; = 0.4 



/imHe£' SN+SW 
ki3Mi 



SN+SW 



Pe=0.118 



,1/2 0I/2 ^1/2 
^ P S.^ SN+SW -"^ SN+SW 

Ri 



and 



: 0.296 



A/lp^l^ d2 
^ ' SN+SW 



(8) 



(9) 



(10) 



The central temperature is thus linearly proportional to the 
ratio of the mechanical energy thermalization efficiency to the 
mass loading factor, and independent of the size of the star- 
burst region 7?* or participation factor Q. Thus it is an excel- 
lent diagnostic tool for quantifying the effects of feedback, 
provided that this temperature can be accurately measured. 

The terminal velocity of the wind fluid 



Voo = 



SN+SW 



\I3M 



SN+SW 



1/2 



(11) 



is also a function of e/(3. An observational determination 
of v'oo is beyond the capabilities of current X-ray observato- 
ries, but the X-ray Microcalorimeter Spectrometer (XMS) in- 
strument on the future International X-ray Observatory (IXO) 
will provide the high resolution non-dispersive X-ray spec- 
troscopy necessary to directly measure the velocity of the hot 
phases in superwinds. 

3.3.4. The metal abundance within the wind fluid 

The mixed SN ejecta and stellar wind material will be en- 
riched in heavy elements compared to the cooler ambient ISM 
within the starburst region. As the emissivity of the thermal 
X-ray emission from a hot plasma depends strongly on the 
metal abundance of the hot plasma it is necessary to consider 
both SN-related enrichment and dilution via mass loading. 

If we assume mass-loading is associated with efficient mix- 
ing, then the elemental abundance Z, wf of an element / in the 
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wind fluid is 

„ Z,-.SN+SW C-^SN+SW + Z,-,coldA^cold Z/.SN+W + (/3 - 1 )Z,xold 

WF = '■ '■ = T. 

C^/,SN+SW + Mcold P 

(12) 

where Z, sn+sw is the abundance of the element in the SN plus 
stellar wind ejecta and Z, coij is the abundance in the ambient 
ISM, and we have made use of Equation|7]which defines 

The primary elemental coolants for plasmas with tem- 
perature T > 10^ K are oxygen, neon and iron (e.g. see 
[Sutherland & Dop"italll993h . Iron is of particular importance 
for this work given the E ~ 6.7 keV and E ^ 6.9 keV lines 
from helium-like and hydrogen-like iron ions. Elements such 
as S, Ar and Ca also produce weaker X-ray lines at energies 
E > 2 keV that may also serve as a diagnostic of the wind- 
fluid 

iLimongi & Chieffil (l2007h present the IMF-averaged SN 
ejecta (both SN II and SN Ib/c) element abundances with re- 
spect to Solar^. Their stellar evolution models include mass 
loss due to stellar winds, but the quoted yields do not include 
the hydrogen and lighter elements returned to the ISM by the 
winds. They find Zpe^sN I'Z.Fe.Q, Zq.sn ~ 6Zo,0, ZNe,SN ~ 

lOZArf.0, Zs^SN ~ 5Zs^0, Zai.SN ~ 4Z^;.^q and Zca.SN ~ SZca,©. 

Note the low O/Fe abundan ce ratio comp a red to the IMF- 
averaged yields presented in iGibson et al.l (11997!) . which is 
most probably a consequenc e of wind mass loss in th e stel- 
lar evolution models used by lLimongi & Chieffil (|2007|) . 

The time-averaged nature of an IMF-averaged yield is ap- 
propriate for continuous star formation with an age of f > 20 
Myr, but the abundances in the ejecta f or a younger or in- 
stantaneous burst may differ from that in iLimongi & Chieffil 
( |2007|) . To investigate this, and the effect of using a different 
set of SN yields, we computed the abundances of the com- 
bined SN ejecta and stellar wind material from the yields 
given by Starburst99. Starburst99 includes wind yields by 
combining the stellar evolutionary wind mass loss rates with 
observed stellar sur f ace ab undances. SN yields are based on 
IWoosley & Weaver! (Il995h for SN type II, but note that the 
the stellar evolution models that lWooslev & Weaverl used did 
not in general include wind mass loss. Starburst99 also as- 
sumes that stars that become Wolf-Rayet stars (initial mass 
M^: > 33 Mq) are assumed to explode as SN type Ib/c and con- 
tribute no metals to the final yields. For a continuous star for- 
mation model the SN plus stellar wind ejecta abundances pre- 
dicted by Starburst99 are relative constant for ages f > 15 Myr 
after the onset of the burst, at levels of Zpe.sN+sw ^ ^^Fc.q, 
Zo,SN+sw ^ 6Zo,0, and Zs.sn+sw ^ 14Zs.q (Starburst99 does 
not calculate the yields of Ne, Ar and Ca). 

Given the uncertainties we will adopt Zpe.sN+sw = SZpe,© 
for these models. For computational convenience we will fur- 
ther assume the same abundance relative to Solar for the other 
elements. 

The stellar and H II region O and Ne abundances i n M82 
are Solar or highe r (fAchterma nn & Lacv 1995; Origli aTt al.l 
12001: ISmith etan i2006b). The are very few measurements 
of iron abun d ances , the only one we are aware of is that of 
lOriglia et al.l (|2004|) who derive a stellar iron abundance of 
Zpe = O.46;!^ j7ZFe.0 from NIR spectroscopy. The total bary- 
onic mass of M82 is ^ 2 x l^O'^M© based on its peak rota- 
tional velocity (see IStrickland et al.ll2004al) . Galaxies of this 
mass have a mean nebular oxygen abundance of ~ 1 .2Zo.0 

^ All abundances quote d in this paper are relative to the scale defined in 
lAnders & Grevessd fl989l) . 



(on the lAnders & Grevessdl 19891 scale), based on the galaxy 
mass-metallicity relationship (iTremonti et al.ir2004 ). Thus the 
measured metallicity of the stellar population and the warm 
ISM in M82 are not unusual given its mass, nor do they sug- 
gest that the ambient ISM has been significantly enriched with 
nucleosynthetic products from the current starburst. We adopt 
•Zpccoid = IZpe,© for our models, again assuming the same 
value for all other elements. 

3.4. Broad-band cooling and X-ray line emission 

The luminosity per unit volume C ( erg s"' cm"^ ) due to the 
thermal emission from an element of gas of density p, tem- 
perature T and metal abundance Z is 

£ = «.«HA(r,Z) = ^"^^ A(T,Z), (13) 

where xh ~ 0.7 1 is the mass fraction of hydrogen for a Solar 
abundance plasma, niu is the mass of a hydrogen atom and 
A(r,Z) is the cooling function for the specific energy band or 
line transition being considered. 

For broad-band X-ray emission, or the total cooling rate, 
the cooling function for a hot optically thin plasma may be 
split into a metallicity-independent component Ah+hcC?^) and 
metallicity-dependent component Az{T) (evaluated for Solar 
abundance), such that A(T,Z) = Au+Hts(T)+Z x AziT), where 
Z is the metal abundance with respect to the Solar This as- 
sumes that the relative metal abundances are Solar, as other- 
wise the contribution from each element must be considered 
separately. 

To assess the total cooling rate w e use the zero-field 
non-equilibrium cooling functions from^ Sutherland & Dopital 
(|T993). We calculated the metallicity-dependent cooling Az 
from their published Solar and zero metallicity cooling func- 
tions. At temperatures above T ~ 5 x 10*' K (typical of the 
vast majority of the models we will consider) the difference 
between equilibrium and non-equilibrium cooling functions 
is negligible. 

To calculate the cooling functions for the E = 2 — 8 keV 
hard X-ray band and the Hea and Lya like lines of S, Ar, 
Ca and Fe we used version 1.3.1 of APEC hot plasma code 
(which assumes collision ionization equilibr ium), accessed 
using the X SPEC spectral fitting program (lArnaudI 119961 : 
ISmith et ani2(j01 ). For the line emission from any specific 
element ; the cooling function is an exactly linear function of 
the metalUcity Z,-, i.e. A(r,Z,) = Z,- x Az(T). 

In Fig. |4] we plot the total and 2-8 keV cooling func- 
tions (split into their metallicity-independent and metallicity- 
dependent components) that we use in our calculations. The 
line-specific cooling functions for hydrogen-like and helium- 
like sulphur, argon, calcium and iron for unit metal abundance 
are also shown. 

The X-ray line emissivities we use apply to hot gas in col- 
lisional ionization equilibrium (CIE). In our the calculations 
using the scaled CC model we evaluate the following expres- 
sion in each radial shell as we numerically integrate outwards: 



If the ionization paramete r ri > 10 '^cm ^s then the plasma is 
Ukely to be in CIE (Mas ai et al.ll2002h . If rj < lO'^cm"^ s 
then the outflow ing gas is ov er-ionized and in recombining 
non-equilibrium. [Masai et all show that the broad band X-ray 
emissivity of such a recombining plasma is less than that of a 
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log T (K) log T (K) 

Fig. 4. — (a) The broad band cooling functions used in the model calculations. The total cooling due to metals for unit metal abundance with respect to Solar 
abundance (1) and cooHng from hydrogen and helium only (2) for non-equilibrium cooling are both fromlSutherland & Dopita" '1993'). The APEC plasma code 
(Smith etal. 2001) was used to calculate the emission in the hard X-ray energy band (£ = 2 - 8 keV) from metals (3) and hydrogen plus helium (4), assuming 
coUisional ionization equilibrium, (b) X-ray cooling curves per unit metal abundance for helium-like and hydrogen-like ions of S, Ar, Ca and Fe in comparison 
to the 2-8 keV cooling curve for all metals (3). 



CIE plasma. The ionization parameter rj scales with thermal- 
ization efficiency and mass loading proportionally as 
as increasing the mass loading or decreasing the thermaliza- 
tion efficiency acts both to increase the plasma density and 
decrease the outflow velocity. 

For parameters appropriate to M82 we find that the plasma 
is typically in CIE at radii r<R^, even when e = /? = 1 . Even if 
we integrate the emission out to radii ^ R^, the total luminos- 
ity is dominated by emission from within justifying our 
use of CIE emissivities. The radial variation in rj and the cu- 
mulative emission integral EI{< R) shown in Fig [3]: demon- 
strate this for one set of possible model parameters. This re- 
sult is consistent with the alternative ionization equilibrium 
and Coulomb relaxation timescale calculations presented in 
Paper 1, which also indicated that the hot plasma within the 
starburst region should be in ionization equilibrium. 

For a single temperature plasma the luminosity is the prod- 
uct of the emission integral EI with the cooling function 
A. The cumulative volume emission integral is defined as 

rR 9 

EI(< R) = Jy 4-Tmi,niiR dR. For a CC-style radial wind this 
is dominated by material within the relatively dense starburst 
region, as shown in Fig[3};- Thus the total emission integral is 
proportional to 

EI oz p V^(x ——-^ , (15) 

en* iiSN+SW 

which again increases as either the mass loading factor /3 is 
increased, or as the thermalization efficiency e is decreased. 
However such variations in /? and/or e are also associated with 
a decrease in the central temperature, and for increased mass 
loading a decrease in the mean metal abundance, so that the 
net luminosity can ultimately decrease even as increasing (3 or 
decreasing e monotonically increase the net emission integral. 

3.5. Observational constraints 

The primary aim of this project is to assess the range of 
starburst properties that are consistent with the hard X-ray Fe- 
line-emitting plasma in M82, and which models are deemed 
acceptable depends on which constraints are used — this is 
unavoidable. These are not an arbitrary set of observational 
and conceptual constraints but have been chosen based on re- 
alistic physical expectations. 

Based on the results and discussion presented in § |2] we 
adopt the following as the full observational constraints on 
the properties of the wind fluid within a radius of 500 pc of 
the center of M82. 



• The total X-ray luminosity of the wind fluid in the E = 2 

- 8 keV energy band must be less than or equal to the 
broad band luminosity estimated from the ACIS-S ob- 
servation for the diffuse hard X-ray emission (Table [T]i, 
i.e. logL2-8keV < 39.65. 

• The Fe Hea line luminosity must be consistent with the 
observed value given in Table [T] In order to be con- 
servative in our estimation of the constraints the ob- 
servations place on the thermalization efficiency and 
mass loading, i.e. to accept the largest number of mod- 
els matching the observational constraints, we adopt 
the weakest constraints on the allowed value of LpeHea, 
i.e. using the ACIS-S observation rather than the bet- 
ter constrained values from the ACIS-I or XMM-Newton 
data sets. Models producing Fe line luminosities within 
68.3% confidence region from the ACIS-S observation, 
37.96 < logLpeHea < 38.39, will be judged as matching 
the constraint on the E ^ 6.7 keV line emission. This 
measurement of LpeHea is slightly larger than, but still 
consistent with the line luminosity estimated from the 
AClS-1 and XMM-Newton observations. 

• The Fe Lya luminosity must be less than or equal to the 
observationally-determined upper limits, of which the 
best determined value is from the XMM-Newton obser- 
vation. We specify that LpeLya /ipeHea should be < 0.34 
for a model to be judged to be consistent with the ob- 
servational constraints (Table [U. 

• The luminosity of the predicted Hea and Lya-like lines 
of S, Ar and Ca must be less than or equal to the ob- 
served values and/or 99% upper limits. Here the best 
observational constraints are from the XMM-Newton 
data (Tabled. Acceptable models must then satisfy all 
of the following: logLsHea < 38.19, logLsLya < 37.43, 

logLArHea < 37.53, logLAiLya < 37.26, logLcaHea < 

37.27, and logLcaLya < 37.14. Upper limits are used 
as the S, Ar and Ca-line-emitting plasma is not neces- 
sarily the same phase as the Fe-line-emitting plasma. 

Models that (a) have total radiative losses less than 30% of 
the net energy injection rates (see AppendixlAl) and (b) predict 
X-ray luminosities consistent with the full set of constraints 
given above are considered as having successfully matching 
the full set of constraints. 
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We will consider the effect of using a different set of ob- 
servational constraints. These alternatives deliberately ignore 
what we know about the iron line emission (the iron lines are 
not used as constraints) and instead attempt to fully account 
for all of either the broad-band hard X-ray luminosity or the 
detected S, Ca and Ar line luminosities: 

• The total X-ray luminosity of the wind fluid in the E = 2 
- 8 keV energy band must be within the 3ct envelope 
given in Table [T] for the model to considered success- 
ful, i.e. 39.59 < logL2-8keV < 39.71. In addition the the 
total radiative losses must be less than 30% of the net 
energy injection rate. We shall refer to these constraints 
as the broad-band-only constraints. 

• The predicted line luminosities for the S Hea, S Lya, 
Ar Hea and Ca Hea must be within ±3cr of the values 
given in Table [3] for the best fit to the XMM-Newton 
spectrum. The allowed limits on the Ar Lya and Ca 
Lya line luminosities are the same as given for the full 
set of constraints. In addition the model must predict 
broad band and total radiative energy losses less than 
the upper limits also used in the full set of constraints: 
logL2-8keV < 39.65 and total radiative losses less than 
30% of the net energy injection rate. We shall refer to 
these constraints as the SArCa constraints. 

We allow larger ranges in the broad-band or line luminosity 
(±3cr) as we found that in practice that it is relatively hard to 
find regions of parameter space where models can satisfy the 
broad-band-only or SArCa constraints. 

4. SIMULATIONS AND THE RESULTS THEREOF 
4.1. Acceptable values of t and [3. 

Mass loading factors of /3 ~ 10 or thermalization efficien- 
cies of e ^ 0.1 mark the practical upper and lower limits on 
these parameters for a M82-like starburst. Each of these val- 
ues on its own will reduce the central temperature to 7^ 10^ 
K, a temperature at or below which the emissivity of the Fe 
Hea emission drops dramatically (Fig.Hji. 

To explore this region of the mass loading and thermaliza- 
tion efficiency parameter space we adopt the following val- 
ues for the other important model parameters: C, = \, a. disk- 
like starburst with c/* = 750 pc and K = 105 pc {R^, = 300 pc), 
ZsN+sw = 5 Zq and Zcoid = 1 Zq. We treat the starburst as con- 
tinuous star formation at a rate of iS = 4.7 Mq yr"' beginning 
30 Myr ago, which leads to -Esn+sw = 3.1 x 10^^ ergs"' and 
A^SN+sw = 1.4M0 yr"'. We shall refer to this combination of 
input model parameters as parameter set A (see Table|5]l. 

We initially computed models covering the full range of 
e : (3 parameter space. To prevent the total number of sim- 
ulations from becoming excessive we considered only a lim- 
ited number of distinct values for e and (3, specifically e = 
and (3 = 2^, where N ranges from to 3 in increments of 
0.25, resulting in the computation of 169 models using the 
2-dimensional hydrodynamical code described in § 13.11 

We find that only limited regions of parameter space yield 
models that satisfy any one of the sets of observational con- 
straints we chose to apply (Fig.|5]l. Furthermore, the regions in 
e:/3 parameter space satisfying the different sets of constraints 
are well separated from one another 

Models that satisfy the full set of observational constraints 
(labeled F in Fig.|5]l have moderately high thermalization effi- 
ciencies (log e > -0.3) and moderately low but non-negligible 



mass loading (0.15 < log/3 < 0.3). These models reproduce 
the observed Fe line fluxes and limits without overproducing 
S, Ar or Ca line emission, but in doing so they can only ac- 
count for a small fraction 10 - 25%) of the total observed 
diffuse luminosity of logLa-skeV ^ 39.6 in the £ = 2 - 8 keV 
energy band (Fig.Et). 

The range of central thermal pressures found in these 
matching models, P^/k = (1.4 - 2.5) x 10^ K cm"-' (given in 
Table|6l) is consistent with estimates of the pressure within the 
starburst region obtained from optical spectroscopy: P/k ^ 
(0.5 - 31 X 10^ K cm"3 (|Heckman et all [19901: ISmifli et alJ 
l2006bHWestmoquette et alJl2007h . 

To match the observed S, Ar and Ca line fluxes using a wind 
model (the region labeled SArCa), models require a lower 
thermalization efficiency and a comparable or slightly higher 
mass loading factor (-0.7 < e < -0.4, log/3 0.4). 

To reproduce the observed diffuse X-ray luminosity in the 
£■ = 2 - 8 keV energy band with thermal emission alone 
requires both high thermalization efficiency and high mass- 
loading factors: loge > -0.3 and (3 >5 (the region labeled 
BB in Fig. |5]l. Models in this region of parameter space fail 
to match the full set of observational constraints because they 
typically predict more S, Ar or Ca emission (the S Lya lumi- 
nosity is shown in Fig. [5]:) than is observed when the Fe line 
luminosity is correctly reproduced, or too little Fe Hea emis- 
sion when they predict less than or equal to the amount of S, 
Ar or Ca line emission observed. 

Thus we can only account for < 50% of the observed E = 2 
- 8 keV diffuse X-ray luminosity with R = 500 pc with the 
sum of the predicted thermal emission from the starburst re- 
gion and the Inverse Compton X-ray emission (Paper I). We 
conclude that either our estimate of the Inverse Compton 
luminosity is in error, or there must be an additional non- 
thermal or quasi-thermal source of har d X-ray continuum 
emission within the starburst region (e.g. iDogiel et alJ[2002l: 
Masai et al. 2002). Solving this enigma will require higher 
resolution spectroscopy to establish the shape of the contin- 
uum and the ionization state of the line emission. 

Models where the total radiative losses are higher than 30% 
of the energy injection rate, and the adiabatic assumption is 
no longer valid, are labeled NA in Fig.|5] Note that for these 
models the temperatures and X-ray luminosities we predict 
are not valid. The location of this region in e : /3 space in 
consistent with the discussion in AppendixlAl For M82's star 
formation rate and starburst region size we find that as long as 
C,0^ /(^ < 1000 then the adiabatic assumption is applicable to 
either a disk-like or spherical starburst. 

Parameters that give rise to highly radiative winds can be 
rejected even though we can not predict their exact X-ray lu- 
minosities. Cooling is so efficient under such conditions that 
the central temperature Tc <C 10^ K and there is no possibility 
of Fe Hea emission (e.g. iTenorio-Tagle et al.ll2007i) . 

To further refine our knowledge of the exact region of the 
allowed e : (3 parameter space for M82's starburst we ran a 
further set of full hydrodynamical simulations a subset of the 
full e : (3 grid described above. We used a linearly-spaced 
grid, varying e between 0.1 and 1.0 in steps of Ae = 0.05 and 
(3 between 1.0 and 3.0 in steps of A/3 = 0.1 (a total of of 494 
models). 

The predicted regions of allowed parameter space for pa- 
rameter set A are shown in Fig. |6l This figure also demon- 
strates that the exact shape of the allowed region using the full 
set of constraints is most-strongly influenced by the combina- 
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Fig. 5. — Predicted luminosities and temperatures as a function of theinput thermalization efficiency e and mass loading factor /3, as predicted using the 
2-dimensional hydrodynamical simulations with parameter set A (see § 14.11 . (a) The logarithm of th e Fe Hea luminosity. The grey-shaded region labeled F 
shows the region of parameter space that matches the full set of observational constraints discussed in § 13.51 The region labeled BB denotes the region satisfying 
the broadband-only constraints, while the region labeled SArCa satisfies the S, Ar, and Ca line constraints. The region labeled NA denotes the region of parameter 
space where the adiabatic assumption ceases to the valid, (b) As (a), except the labeled contours show the logarithm of the luminosity in the £ = 2 - 8 keV energy 
band, (c) As (a), except the labeled contours show the logarithm of the S Lya line luminosity, (d) As (a), except the contours show the logarithm of the central 
plasma temperature Tc in Kelvin. 



tion of the constraints on the Fe Hea line and the upper limit 
on the S Lya line**. For a given value of the thermalization ef- 
ficiency e the lower boundary of acceptable mass loading fac- 
tors /? is set by requirement to achieve sufficient Fe Hea emis- 
sion (the necessary condition is 37.96 < logLpeHea < 38.39), 
while the maximum /3 is set by the requirement not to exceed 
the observed S Lya luminosity (logLsLya < 37.43). 

For a given value of e the highest allowed fi is a less likely 
solution than a lower allowed value of /3. The highest allowed 
values of /3 produce the entire observed S Lya flux in the wind 
fluid, without leaving leeway for S Lya emission from the soft 
X-ray-emitting plasma within the starburst region. 

Were we to change the S Lya constraint from an upper limit 
to a specific range in values around the observed line luminos- 
ity the resulting allowed region would be a narrow strip along 
what is currently the upper fi margin of the full allowed re- 
gion. Adopting such a constraint would be somewhat ad-hoc, 
given that the S Hea and Ar Hea lines are equally well de- 

* For this reason we have not provided plots of the S Hea, Ar or Ca line 
luminosities as a function of e and /3. 



tected as the S Lya line (Table [3]). Furthermore, an attempt 
to match the S Hea, Ly a and Ar Hea line fluxes, rather than 
treating them as upper limits, simultaneously with the Fe Hea 
line flux leads to no acceptable models being found. Ulti- 
mately this is because the S, Ar and Ca line fluxes and upper 
limits are indicative of a plasma temperature logT < 7.2 (see 
Table|2]or the SArCa region on Fig.|6}l), yet such a low tem- 
perature plasma produces too little Fe Hea emission to match 
the observations. 

As shown in Fig. |7] these same general regions of accept- 
able parameter space are reproduced in simulations making 
use of the scaled CC model, although in detail there are mi- 
nor difference caused by the < 30% errors in the luminos- 
ity predicted by the scaled CC model compared to the multi- 
dimensional hydrodynamical simulations (AppendixlBl). 

4.2. Other input model parameters. 

Given the computational cost of the full hydrodynamical 
simulations it is preferable to use the scaled CC model to ex- 
plore the effect of varying other model parameters such as 
starburst region size, metal abundance and participation fac- 
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TABLE 5 

Combinations of starburst model parameters explored in the 
simulations. 



Parameter Set 


SF Mode/Age 


£sN+SW 


^SN+SW 


R* 


2sN+SW 


c 






( erg s"' ) 


(M0yr-l) 


(pc) 


(Z©) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


U) 


A 


CSF / 30 Myr 


3.1 X 10*2 


1.4 


300.0 


5.0 


1.00 


B 


CSF/ lOMyr 


1.9x 10« 


1.1 


300.0 


5.0 


1.00 


C 


SIB / 10 Myr 


4.5 X lO''^ 


2.2 


300.0 


5.0 


1.00 


D 


CSF / 30 Myr 


3.1 X 10*^ 


1.4 


125.5 


5.0 


1.00 


E 


CSF / 30 Myr 


3.1 X ltf*2 


1.4 


404.9 


5.0 


1.00 


F 


CSF / 30 Myr 


3.1 X W*^ 


1.4 


300.0 


2.5 


1.00 


G 


CSF / 30 Myr 


3.1 X 10*^ 


1.4 


300.0 


7.5 


1.00 


H 


CSF / 30 Myr 


3.1 X ltf*2 


1.4 


300.0 


5.0 


0.71 


I 


CSF / 30 Myr 


3.1 X W*^ 


1.4 


300.0 


5.0 


0.50 



Note. — Column 1 : Adopted name of combinations of model parameters. Column 2: The mode of 
star formation assumed, either continuous star formation (CSF) or an single instantaneous burst (SIB) 
beginning or occurring the specified time before the present. Columns 3 and 4: The mechanical energy 
and mass injection rates (See § |3.3.H associated with the assumed star burst. Column 5: The effective 
starburst radius adopted in calculations using scaled CC model. Column 6: The assumed metal abun- 
dance of heavy elements relative to Solar (see § |3.3.4t . Column 7: The assumed participation fraction 
(see § |3.3.3| or Appendix[cj. 




Fig. 6. — As Fig. [5] except using a finer grid of models distributed linearly in a smaller region of e : /3 parameter space. The numbered contour levels are of log 
LfcHcc (a), log L2-8keV (b), log LsLya (c) and log Tc (d). 

tor. Our experience with parameter set A indicates that the to an accuracy of ^0.1 for the full set of observational con- 
scaled CC model reproduces the allowed values of e and [3 straints that include the Fe lines. 
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Fig. 7. — (a) A comparison between the region of allowed e : /3 parameter space for M82 predicted by the scaled CC model (shaded regions) in comparison to 
the regions predicted by the full two-dimensional hydrodynamical simulations (regions enclosed by dot-dashed lines, previously shown in Fig. [5), using parameter 
set A. (b) As in (a), except the comparison uses the finer giid of models distributed linearly in a smaller region of e : /3 parameter space. 
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Fig. 8. — (a) A comparison between the region of allowed e : (3 parameter space for M82 predicted by the scaled CC model using parameter set B (shaded 
regions) in comparison to the equivalent regions for parameter set A (dashed lines), (b) As in (a), except the comparison is between parameter set C and parameter 
set A. Parameter sets B and C differ from A in the assumed starburst star formation mode and age. 



Parameter sets B through I are shown in Table |5] These al- 
ternate parameter sets explore the influence of the chosen star 
formation mode and age, starburst region size, supernova and 
stellar wind metal abundance, and participation factor on the 
region of allowed e : (3 parameter space. All of these calcula- 
tions were performed using the scaled CC model alone. Pre- 
dicted central density, pressure and temperature, along with 
the terminal velocity, total cooling, broad-band hard X-ray 
(£ = 2 - 8 keV band) luminosity and Fe Hea and S Lya line 
luminosities for e = 1 and /3 = 1 are shown in Table |7] 

Parameter set B differs from the previously discussed pa- 
rameter set A only in that is assumes that ongoing star for- 
mation began 10 Myr ago, rather than 30 Myr ago. As the 
younger stellar population is intrinsically more luminous, less 
star formation is required to match the IR luminosity of M82 
and consequently the mechanical energy and mass injection 
rates in parameter set B are lower than in parameter set A. 
Furthermore the ratio of energy to mass injection is lower, re- 
sulting in a lower central temperature of Tc ~ 8 x 10^ K for 
unit e and /3 compared toT^^ 10** K in parameter set A, and a 
lower temperature than parameter set A for any given value of 



e and f3. Cumulatively these differences result in a significant 
shift in the location of the the models in e : /? parameter space 
that match either the full or SArCa constraints (see Fig.[8t). 

Generation of plots equivalent to Fig.|6](not shown) reveals 
that the region matching the full set of constraints is again 
determined by the intersection of the Fe Hea and S Lya ob- 
servational constraints. Furthermore, at a given value of e and 
/3 the predicted broad band and (in particular) line luminosi- 
ties are lower for parameter set B than for parameter set A, 
with the least difference in luminosity occurring at e = /3 = 1 . 

Application of Equation [15] indicates that parameter set B 
yields emission integrals consistently ^0.1 dex lower than 
parameter set A for all values of e and f3. For the Fe and S 
lines the lower temperature of parameter set B when (3 
results in a higher cooling function that offsets the reduced 
emission integral because the lower temperature is closer to 
the peak of the line-specific cooling functions, such that the 
predicted values of LpeHea and LsLya are comparable to those 
of parameter set A (Table |7]). But over much of the e : [3 pa- 
rameter space the lower temperatures in parameter set B re- 
sult in lower line emissivities that in combination with the 
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Fig. 9. — (a) A comparison between the region of allowed c : fi parameter space for M82 predicted by the scaled CC model using parameter set D (shaded 
regions) in comparison to the equivalent regions for parameter set A (dashed lines), (b) As in (a), except the comparison is between parameter set E and parameter 
set A. Parameter sets D and E differ from A in adopting respectively smaller and larger starburst effective radii R*. 





Fig. 10. — (a) A comparison between the region of allowed e : j3 parameter space for M82 predicted by the scaled CC model using parameter set F (shaded 
regions) in comparison to the equivalent regions for parameter set A (dashed lines), (b) As in (a), except the comparison is between parameter set G and parameter 
set A. Parameter sets F and G differ from A in adopting respectively smaller and larger supernova and stellar wind metal abundances Zjn+sw- 
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Fig. 1 1 . — (a) A comparison between the region of allowed e : /3 parameter space for M82 predicted by the scaled CC model using parameter set H (shaded 
regions) in comparison to the equivalent regions for parameter set A (dashed Unes). (b) As in (a), except the comparison is between parameter set I and parameter 
set A. Parameter sets H and I differ from A in adopting lower values of the participation factor f = 0.71 and C = 0.50 respectively, i.e. effectively lower star 
formation rates. 
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TABLE 6 

SUPERWIND PROPERTIES THAT MATCH THE FULL SET OF OBSERVATIONAL CONSTRAINTS. 



Parameter Set Wmatch Range in parameter consistent with full constraints 



(1) (2) (3) (4) (5) (6) (7) (8) 



A" 52 0.55- 1.0 1.4-2.2 7.56-7.82 7.15-7.38 34.31 -34.54 1555-2097 

A 75 0.40- 1.0 1.2-2.3 7.51 -7.84 7.08-7.40 34.23-34.55 1470-2164 

B 26 0.75- 1.0 1.8-2.3 7.49-7.66 7.13-7.24 34.29-34.39 1434- 1745 

C 32 0.45- 1.0 1.0- 1.4 7.64-7.87 7.22-7.47 34.37-34.62 1709-2239 

D 21 0.45- 1.0 1.0- 1.4 7.67-7.87 7.80-8.05 34.20-34.44 1778-2240 

E 60 0.55- 1.0 1.7-2.6 7.48-7.76 6.94-7.17 34.35-34.58 1426- 1976 

F 47 0.60- 1.0 1.9-2.6 7.50-7.72 7.26-7.43 34.41 -34.58 1452- 1874 

G 76 0.30- 1.0 1.0-2.0 7.50-7.87 6.96-7.37 34.11 -34.52 1452-2240 

H 43 0.65- 1.0 2.1 -2.8 7.47-7.70 7.15-7.29 34.30-34.44 1409- 1829 

I 



Note. — Column 1: Parameter set name. Column 2: Nmmch is the number of models that match the full set of observational 



constraints, out of the 494 models calculated for each parameter set. Columns 3 and 4: The range of thermalization efficiencies 
€ and mass loading factors /3 found to match the full set of observational constraints. The remaining columns show the range in 
a variety of physical properties of the M82 superwind found in the models that match the full set of observational constraints: 
central temperature Tc (column 5, in Kelvin), central pressure P^^/k (column 6, in units of K cm~^), rate of change of superwind 
momentum /j^v'f ~ ''oo (column 7, units g cm s~") and wind terminal velocity v^d (column 8, units km s~' ). 

^ Calculated using 2-dimensional hydrodynamical simulations. All other parameter sets were calculated using the scaled 1- 
dimensional.Chevalier & Clegg. model. 
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lower emission integral result in lower line luminosities than 
in parameter set A. Ultimately the regions in e : /? parameter 
space that match the full and SArCa constraints in parame- 
ter set B require higher mass loading and higher thermaliza- 
tion efficiency than in parameter set A in order to increase the 
emission integral without further decreasing the temperature. 
Increasing (3 strongly increases the net emission integral as 
EI (X (3^ /e, although it slightly decreases the plasma temper- 
ature (which must then be offset by increasing the thermal- 
ization efficiency). In contrast, reducing the thermalization 
efficiency is counter productive as the increase in EI is coun- 
teracted by the decrease in temperature that consequently re- 
duces the line emissivities. 

Parameter set C represents the star formation in the nucleus 
M82 as a single instantaneous burst occurring 10 Myr ago. 
Compared to parameter set A this results in a higher mechan- 
ical energy and mass injection rate, although the ratio of en- 
ergy to mass injection is much closer to that in parameter set 
A than the lower values of parameter set B. 

The higher mass injection for parameter set C leads to 
higher predicted X-ray luminosities for a given value of e and 
(3 than for parameter set A, which has a significant effect on 
the regions of e : (3 parameter space that match the various 
constraints (Fig.lHJi). The region matching the E = 2-2, keV 
broad band luminosity appears at much lower f3 than in pa- 
rameter set A, although in terms of absolute mass injection 
Mtot the location of this region is consistent between parame- 
ter set A and parameter set C. 

Similarly, lower values of /3 are required to match the 
full set of constraints, the allowed region extending between 
0.4 < e < 1 and 1 < /3 < 1 .4. Models with e = /3 = 1 are again 
excluded as possible matches to the observed properties of 
M82, but for different reasons than in the case of parameter 
sets A and B. For parameter set C models with f3 ^ 1 and 
e ^ 1 generate sufficient Fe Hea luminosity to match the re- 
quired constraint, but these models predict a Fe Lya/Hea line 
ratio greater than the upper limit of 0.34 (i.e. Tc is too high), 
and hence do not match the full set of constraints. A combi- 
nation of constraints shapes upper edge of the allowed region. 
In addition to the requirement not to overproduce the S Lya 
emission that we encountered with parameter sets A and B, 
the upper edge of the allowed region in parameter set C is 
also set where the Fe Hea luminosity exceeds the observa- 
tional limits. 

Parameter set D adopts the compact starburst region dis- 
cussed in § 13.3.21 In terms of the effective radius of the 
starburst region parameter set D uses = 125.5 pc instead 
of Ri, = 300 pc, but is otherwise identical to parameter set 
A. Consequently the gas density within the starburst region 
(hc = 0.91 cm"^ ) is considerably higher than in parameter set 
A («c = 0. 16 cm"-* ) or even parameter set C (Wc = 0.26 cm"^ ), 
and results in a significant increase in the X-ray emission de- 
spite the reduction in starburst region volume (see Equ. [TsT l. 
The regions of e : (3 parameter space over which parameter 
set D matches the various sets of constraints is very similar 
to that parameter set C (Fig. |9^) for the same reasons we 
discussed above. The region of e : /3 parameter space that 
matches the full set of constraints is smaller than in parame- 
ter set C because of the greater increase in luminosity caused 
by the higher gas density, and no region of e : /3 matches the 
SArCa constraints. It is possible that with tighter observa- 
tional constraints on LpeHea and LsLya there would be no re- 
gion of e: /3 parameter space that would match the full set of 



observational constraints. 

The lower volume of the starburst region adopted in pa- 
rameter set D also leads to higher central pressures. Mod- 
els that match the full set of observational constraints have 
7.80 < logPjk < 8.05, values that are 2 - 4 times higher 
than the highest independent observational estimates of the 
thermal pressure in the core of M82. This leads us to suspect 
that parameter set D is not a valid physical representation on 
the M82 starburst. 

Parameter set E adopts the largest of the possible starburst 
region geometries discussed in § 13.3.21 for which the effective 
radius is = 406 pc and the resulting central number density 
is He = 0.09 cm""* . The region that matches the full set of ob- 
servational constraints now moves toward higher e and higher 
(3 with respect to the same region in parameter set A (Fig.|9j3), 
for similar reasons to those we discussed regarding parameter 
set B. Increased mass loading is required to make up for the 
reduction in gas density due to the larger starburst, but the 
thermalization efficiency must also increase so as to maintain 
a central temperature that maximizes the Fe Hea emissivity 
(7.5 < logTc < 7.8, for parameter sets A, B and E. See Ta- 
ble©. 

In § 13.3.41 we discussed the expected metal abundance of 
the merged SN ejecta and st e llar w ind ejecta. The calcu- 
lations of Limongi & Ch ieffj (|2007|) predict a higher iron 
abundance for the SN ejecta (Zpe sn = 7ZFe.(T)) than the Star- 
burst99 calculations that combine IWoosley & Weavedl (Il995h 
SN yields wi th stellar wind ejec ta (Zpe.sN+sw = 4ZFe.0). Fur- 
thermore the iLimongi & Chieffil calculations predict that sul- 
phur, argon and calcium are less enriched than iron, while 
Starburst99 predicts that sulphur is more strongly enriched 
than iron. In all of the calculations we have discussed until 
now we adopted a uniform relative enrichment of all elements, 
ZsN+sw = 5Zq . Parameter sets F and G explore the influence 
of the adopted metal abundance on the region of e : /3 param- 
eter space that match the observational constraints: parameter 
set F adopts Zsn+sw = 2.5Zq while parameter set G adopts 
ZsN+sw = 7.5Z0. Given the already large uncertainties in the 
yields we feel that exploring the effect of adopting different 
yields for different elements is beyond the scope of this paper 

The lower metal abundance in parameter set F leads to 
lower lower S, Ar, Ca and Fe line luminosities than for param- 
eter set A. For (3=1 the line luminosities are exactly half those 
in parameter set A, but the fractional difference becomes less 
as /? increases and the influence of the mass-loaded material 
on the net metal abundance of the wind fluid Zwf becomes 
larger (Equ. fT2l i. The result in terms of the region of e : /3 
parameter space that matches the observational constraints is 
shown in Fig.fTOh. and is very similar to that in parameter sets 
B and E that also generate lower line luminosities for different 
reasons. 

Conversely the higher metal abundance in parameter set 
G leads to higher line luminosities than parameter set A, al- 
though they are generally less than those in parameter sets C 
and D (see e.g. Table|7|. The e : /3 region of matches is shown 
in Fig.fTOb. The upper f3 edge of the region that matches the 
full set of observational constraints is set by the upper limit on 
the S Lya luminosity, while the lower (3 edge of this region is 
set by a combination of the Fe Lya/Hea line ratio constraint 
and the minimum acceptable Fe Hea line luminosity. 

Parameter sets H and 1 explore the effect of varying the par- 
ticipation fraction C (see § 13.3.3b . Reducing ( has the effect 
of reducing the effective star formation rate without altering 
the energy per particle, and the emission integral and lumi- 
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nosities for a given value of e and /3 are equal to those from 
parameter set A except multiplied by a factor C,'^. For param- 
eter set H, where Q = 2"'/^, the factor 2 reduction in all lu- 
minosities moves the regions that match both the SArCa and 
full set of constraints to higher /3 and higher e (Fig.fTTh). Pa- 
rameter set 1 has ( = 0.5 and the reduction in luminosity is so 
significant that the region matching the SArCa constraints has 
moved off of the region shown in Fig. [TTb to 0.5 < e < 0.9, 
4-5 ^ f3 ^ 5.5. More significantly, there are no values of e or 
P that can match the full set of observational constraints when 
using parameter set 1. 

5. DISCUSSION 

Inl Strickland & HeckmanI (l2007h we concluded that the dif- 
fuse Fe Hea emission from the central ^ 500 pc of M82 could 
have been produced by thermalized SN ejecta and stellar wind 
material from the starburst based on relatively simple esti- 
mates of the plasma properties of the emitting material. 

The more rigorous approach to the interpretation of the dif- 
fuse hard X-ray emission taken in this paper strengthens that 
conclusion: the plasma responsible for the hard X-ray iron- 
line emission is the same gas that drives the larger-scale su- 
perwind in this starburst galaxy. Furthermore our method al- 
lows us to make quantitative statements regarding supernova 
feedback and the properties of M82's starburst-driven super- 
wind. 

5.1. Supernova Thermalization (Heating) Efficiencies and 
Mass Loading 

The primary aim of this paper was to establish what con- 
straints could be made on the efficiency of supernova heating 
and the net rate of mass injection within a starburst region. 

It is generally accepted that in more normal, or quiescent, 
star forming galaxies such as the Milky Way on average only 
^ 10% of the 10^' erg of initial kinetic energy per supernova 
ends up heating or moving th e ISM, with the remainder be- 
ing lo st as radiation (see e.g. Thornton et al.llT998t lEfstathioul 
I2OOOI) . Although it had long been suspected that the efficiency 
of supernova heating is a function of the local interstellar and 
star forming conditions there was no direct method by which 
to measure this thermalization efficiency in a starburst galaxy 
until, as described in § [T] the launch of the Chandra X-ray 
Observatory. 

Theoretical arguments have been made for either lower or 
higher SN heating efficiencies in more actively star-forming 
galaxies. iSteinmetz & Navarro! (Il999h argued that radiative 
energy losses must always be severe becaus e massive stars 
form in regions of high gas density. In contrast lHeckman et al.l 
argue that thermalization must be efficient in star- 
bursts because a high supernova rate per unit volume leads 
to an increased filling factor of the hot low density ISM, 
which in turn reduces the radiative energy losses experi- 
enced by subsequent SNe. Thus the SN thermalization effi- 
ciency depends on the SN-modulated porosity of the ISM (e.g. 
McKee & Ostriker 1977 ). Rece nt numerical simulations by 
Melioli & de Gouveia Dal Pino (2004) appear to qualitatively 
demonstrate such a time-dependent non-linear thermalization 
efficiency. Intriguingly they find that the thermalization ef- 
ficiency (which they term the heating efficiency) can only be 
low (< 10%) or very high (^ 100%) for any significant period 
of time, transitioning between to two states very rapidly. 

Let us briefly review how the hard X-ray line emission 
from a starburst region is related the average thermalization 



efficiency e and the mass loading factor (3. Both of these 
parameters strongly affect the net volume emission integral 
EI cx (Equation [TsTi and the line emissivities A(T,Z) 
through the plasma temperature T (x e/f3 (Equation|8] Fig.|4|i 
and the metal abundance of the wind fluid Zwf (EquationfT2ll. 

In the specific case of M82 it is the Fe Hea and S Lya 
lines that drive the allowed values of e and f3 we derive. For 
many of the sets of model parameters we have explored it 
should be noted that e = /3 = 1 will not reproduce the observed 
constraints because too little Fe Hea emission is produced 
as the plasma density within the starburst region is too low. 
Thus successful models must have a higher volume emission 
integral which implies /9> 1 and/or e < 1 . 

As we described in § 14.11 the lower boundary of acceptable 
mass loading factors /3 for a given value of the thermalization 
efficiency e is set by requirement to achieve sufficient Fe Hea 
emission (37.96 < logLpeHea < 38.39), while the maximum 
/3 is set by the requirement not to exceed the observed S Lya 
luminosity (logLsLya < 37.43). 

Furthermore the central temperature (and hence e/(3) must 
be high enough that the Fe Hea emissivity is significant with 
respect to the emissivity of the lower ionization lines, in par- 
ticular the sulphur lines. This constraint ultimately controls 
the minimum allowed value of the thermalization efficiency 
e, and is the reason why the different parameter sets tend to 
have similar minimum values of T^, specifically logT^ > 7.5 
(see Table. |6l). 

The central temperature can not exceed 7.8 x 10^ K 
(logTc = 7.9) as otherwise the upper limit on the Fe Lya/Hea 
line ratio would be violated. For a given value of e this con- 
straint determines the lowest allowed values of /? in cases 
where the Fe Hea luminosity is sufficiently high enough to 
match the observational constraints. This upper temperature 
constraint is only significant when the model parameters gen- 
erate higher central densities or luminosities (more mass in- 
jection, a more compact starburst, or a higher metal abun- 
dance in the ejecta), and thus there is a broader range in the 
maximum allowed than in the minimum allowed T^. 

Returning to the consideration of our results, we find that 
for all parameter sets that produce at least one valid solu- 
tion that can match the full set of observational constraints 
for M82, those valid solutions always include thermalization 
efficiencies of 100%. 

Furthermore e is never less that 30%, and over the range 
of parameter sets we have explored the minimum valid ther- 
malization efficiency e^in ranges from 30% to 75% (see Ta- 
ble |6]l. Within this range there is a small degree of variation 
associated with different starburst model parameters we have 
explored. We find that the starburst history and ejecta metal 
abundance Zsn+sw have some effect on the allowed thermal- 
ization efficiency, while the starburst region size appears to 
have little influence on the minimum allowed thermalization 
efficiency (Table|6l). For example in a younger continuous star 
formation (CSF) model for the burst (parameter set B) the 
minimum allowed thermalization efficiency is e^in = 0.75 in 
comparison to the older CSF we adopt as the default starburst 
history (e^in = 0.4) or the single instantaneous burst model 
(parameter set C, e,nin = 0.45). 

That £ must be > 30% is a striking result. If the diffuse E ~ 
6.7 keV iron line in M82 arises in a collisional plasma then the 
average thermalization efficiency in the M82 starburst region 
can be as high as 100%, while it can not be as low as 10%. Our 
re sults give quant i tative sup port to the arguments presented 
by lHeckman et alJ (1 19901) and lMeUoU & de Gouveia Dal Pind 
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TABLE 7 

STARBURST region properties for e = 1 AND /3 = 1 . 



Parameter Set 


log«c 


logPc/*: 




Voo 

(kms"') 


rn 
(pc) 


logLtot 


logL2-8ki;V 




logLsLya 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


A 


-0.80 


7.22 


8.02 


2651 


278 


38.67 


38.27 


37.46 


36.22 


B 


-0.85 


7.06 


7.91 


2341 


278 


38.56 


38.17 


37.45 


36.24 


C 


-0.59 


7.40 


7.98 


2548 


285 


39.09 


38.70 


37.92 


36.69 


D 


-0.04 


7.98 


8.02 


2651 


121 


39.17 


38.75 


37.93 


36.74 


E 


-1.06 


6.96 


8.02 


2651 


369 


38.49 


38.09 


37.27 


35.98 


F 


-0.80 


7.22 


8.02 


2651 


278 


38.51 


38.11 


37.16 


35.92 


G 


-0.80 


7.22 


8.02 


2651 


278 


38.80 


38.39 


37.64 


36.40 


H 


-0.95 


7.07 


8.02 


2651 


272 


38.37 


37.97 


37.16 


35.92 


I 


-1.10 


6.92 


8.02 


2651 


266 


38.07 


37.67 


36.86 


35.62 



Note. — Column I : parameter set name. Columns 2, 3 and 4: Logarithms of the central number density ( cm"^ ), the central 
pressure divided by Boltzmann's constant, and the central temperature (K) respectively for models where e = I and /3 = I . Column 
5: Wind terminal velocity. Column 6: The radius at which the ionization parameter 7] (see Equation I 141 drops below lO'^ cm"^ s. 
Columns 7, 8, 9 and 10: Logarithms of the plasma luminosity (ergs^^ ) within a radius of 500 pc. Columns are the total cooling 
luminosity, the £ = 2 - 8 keV energy band X-ray luminosity, the Fe Hea line luminosity and the S Lya line luminosity respectively. 



( |2004|) for a higher thermalization efficiency in starbursts than 
in normal star-forming galaxies. The issue of whether the 
thermalization efficiency must be 100%, or whether high but 
intermediate values are allowed, remains an interesting but 
open question. 

With respect to the mass-loading of the superwind in M82 
we find that the majority of the parameter sets we have consid- 
ered require some centralized mass loading of the wind fluid, 
although at low- to-moderate levels. Over all parameter sets 
1 < /3 < 2.8, but only parameter sets C, D and G allow /? = 1 
and more typically 1-5 < /3 < 2.5. 

Parameter sets C, D and G allow zero additional mass- 
loading (/3 = 1) as a valid solution when coupled with lower 
values of the thermalization efficiency (e 30 - 70%). In pa- 
rameter set C this is because the single instantaneous starburst 
injects more mass than the constant star formation model used 
in the other parameter sets, and hence does not require addi- 
tional sources of mass to increase the plasma density within 
the starburst region. In parameter sets D and G the non-mass- 
loaded plasma within the starburst region is sufficiently lu- 
minous without additional mass-loading because of either the 
compactness of the starburst (D), or the higher metal abun- 
dance (G), increases the net luminosity. High values of e are 
excluded as valid solutions when (3=1 for these three param- 
eter sets because the high central plasma temperature violates 
the Fe Lya/Hea line ratio constraints. 

Thus these results suggest that it is plausible that the raw 
SN and stellar wind ejecta from the massive stars in M82's 
starburst have been efficiently mixed with a roughly equiva- 
lent mass of initially cooler ambient gas, although no mass 
loading is also a possible solution. The net effect is an av- 
erage mass outflow rate from the starburst region in the very 
hot phase of the wind fluid of Mtot ^ 2.5 Mq yr"', approxi- 
mately half of the gas consumption due to the star formation 
rate alone (§ 13.3.11 ). The allowed range of wind fluid outflow 
rate covering all parameter sets except parameter set D lies in 
the range 1.4 < MtatiMe yr"') < 3.6 (a fraction - 0.3 - 0.8 
of the star formation rate for the Salpeter IMF we adopt). 

This level of mass-loading is significantly lower than the 
values others have derived from soft X-ray observations of 
superwinds. Suchk ov" et al.l (|l996) used fairly sophisticated 
1 -dimensional hydrodynamical modeling to estimate that a 
mass-loading factor of /3 ^ 3 - 6 was required in to match the 



soft X-ray properties of M82 known at that time from Einstein 
and ROSAT HRJ observations. For the dwarf starburst galaxy 
NGC 1569 IMartin et alj (|2002|) estimated /3 - (8 - 35) x ry'/^ 
(where r] is the volume filUng factor of the X-ray-emitting 
plasma) based on the estimated ma ss of soft X-r ay-emitting 
gas seen in Chandra observations. iGrimes et al.! (2007) de- 
rived a value of /3 ^ 5 ry'^^ for the soft X-ray-emitting plasma 
in the local Lyr nan Break Ga laxy analog Haro 1 1 . Both 
iMartin et alj and iGrimes et alj obtain similar estimates of P 
based on the best-fit metal abundances in the soft X-ray- 
emitting plasma compared to the expected abundances of a 
SN-ejecta-enriched plasma (see Equation[T2l). 

However a direct comparison between our findings and the 
existing results in the literature is somewhat misleading. 

The results from the literature are all based on soft X-ray 
observations of starbursts that probe gas in the temperature 
range 2-10 million K, not the hotter 30 - 70 million gas 
probed by the observations of M82's diffuse hard X-ray emis- 
sion. The high mass loading factors simply reflect the fact that 
there is more mass in the soft X-ray emitting plasma than can 
be accounted for by SN and stellar wind ejecta alone. 

Soft X-ray mass loading factors are further based on mea- 
surements that cover the entire soft X-ray nebula associated 
with the target starburst galaxy (over ^ 5 kpc-scales for M82 
and Haro 11, and over ^ 1 kpc in NGC 1569), whereas our 
measurement constrains the mass loading of the wind fluid 
within the ~ 500 pc starburst region itself (in the terminology 
introduced by ^S uchkov et al.l (IT996) this is central mass load- 
ing). It is plausible that the wind fluid experiences additional 
mass-loading as it flows outward beyond the starburst region 
(distributed mass loading). 

The older measurements referred to above can not differen- 
tiate between a localized or a global mass-loading process (by 
global we mean a process that is essentially volume filling), 
although we have previously argued that a significant frac- 
tion of the soft X-ray emission must ari se from a low-volume 
filUng factor component of the wind dStrickland et al.ll2000l 
I2002i l2004al) . In contrast the mass loading factor derived in 
this paper is for a volume-filling mass loading. 

It is worth noting that it is not the case that the hottest 
plasma in a superwind is necessarily also the least mass- 
loaded phase. The technique we have presented could detect 
a mass loading factor of ~ 5 - 10 for the very hot gas in the 
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starburst region. The line luminosities are strong functions 
of the mass loading factor (Fig. |5]l, so that both the Fe and 
S line emission would be more intense than observed in the 
case of a higher mass-loading factor Both low thermalization 
efficiency and high mass loading are more readily observable. 
In a sense it is the relative faintness of these lines from M82 
that force us toward the allowed values of the thermalization 
efficiency and mass loading factor given in Table |6] Even 
if part of the starburst region experienced higher mass load- 
ing than another part our final estimate of the allowed mass 
loading factor would be skewed toward the higher value as 
/? is essentially luminosity-weighted. Thus we conclude that 
the T >3 X 10^ K component of M82's superwind (the wind 
fluid) genuinely has a relatively low degree of mass-loading 
within the starburst region. 

5.2. Physical Properties of the M82 Starburst 

Despite the broad range of input model parameters we have 
explored in the various parameter sets we find that the physi- 
cal properties of the base of the M82 superwind in the allowed 
models are all relatively similar (Table|6l). 

We have akeady described the relatively narrow allowed 
range for the central temperature, 7.5 < logTc < 7.9, and the 
physical constraints that drive this result. The result is all the 
more notable given the relatively poor constraints provided by 
the traditional line -ratio diagnostics shown in Table |2l 

Similarly the central pressure and wind momentum injec- 
tion rate are constrained to lie within a range of ^ 0.5 dex 
(7.0 < \ogPc/k < 7.5, 34. 1 < log /^wf < 34.6), if we exclude 
parameter set D from consideration because of its unrealisti- 
cally high pressure. These pressures are in good agreement 
with independent estimates of the thermal pressure within the 
ISM in the starburst region, as previously discussed. 

The momentum injection rate is an important parameter 
with regard to the acceleration of the WNM and WIM by the 
wind fluid (see e.g. Heckmanetal. 1990, 2000)^. By way 
of comparison the momentum injection rate assuming 100% 
thermalization efficiency and no mass loading is log pwF = 
34.2 for the star formation models adopted in parameter sets 
A and B, and logp'wF = 34.6 for parameter set C. Thus taking 
account of imperfect thermalization and additional (central) 
mass loading does not significantly alter the M82 starburst 
momentum injection rate from the values one would estimate 
from SNe and stellar winds alone. 

5.3. Characteristic Velocities of the M82 Superwind 

The moderately tight limits on the central temperature we 
have derived also make for tight limits on the terminal veloc- 
ity of the wind fluid, as oc T^^^ (Equation [TT]). Over all 
the parameter sets, models that match the full set of the ob- 
servational constraints predict Voo in the range 1410 - 2240 
kms"' . 

This estimate of Voo is significantly higher than previous 
estimates of the terminal velocities in local starbursts (includ- 
ing M82) based on the temperature of the soft X-ray emit- 
ting gas, Voo^sx = 440 x (kTsx/0.25keVy/^kms-^ (see e.g. 

' The terminal velocity of a cloud of radius and mass Scloud ™d M^i^ud^ 
initially exposed to a radially-expanding SN-driven wind at a radius Si from 

the origin, is proportional to (/iwp/^i X ^cloud/^^clond^'^"- Notably the termi- 
nal velocity of the cloud can be lower than the terminal velocity of the wind. 
For a radiation-pressure driven cloud the final result is essentially the same, 
with the substitution of p^.^^ = Lbol/c for pwF- Clouds of entrained mate- 
rial are effectively momentum-driven irrespective of whether they are in an 
energy-driven hot wind or a radiation-driven wind. 



lMartinll999l:lHeckman et al.l2000l:lMartinl2005h . However, it 
has been long understood that these estimates were effectively 
lower limits on given the technical difficulties associated 
with robustly detecting a hotter and more tenuous (and thus 
intrinsically fainter) component of the wind. 

Strictly speaking the terminal velocity is the characteristic 
velocity the wind fluid would asymptote to once outside the 
starburst region, if it expanded into a vacuum (see Fig. [3]i. 
For a superwind expanding into a gaseous medium there are 
additional characteris tic velocit ies, some of which are phase- 
dependent ( Stricklan d & Stevens 2000, Strickland & Dinge, 
in preparation), but all are less than or equivalent to Voa- Nev- 
ertheless v'oo remains an accurate measure of both the effec- 
tive energy per particle of the most metal-enriched gas and 
the maximum velocity attainable by gas in a pressure-driven 
superwind. 

For a wind that has broken-out of both the gaseous disk 
and halo of its host galaxy a comparison between Voo and 
the local escape velocity from the gravitational potential Vesc 
should give a qualitative indication of whether the material in 
the wind can escape. 

Estimates of M82's escape velocity (and mass) have sur- 
prisingly large uncertainties for such a well-studied galaxy, 
because of its very unusual rotation curve and the possibility 
that its encounter with M8 1 has stripped it of its dark mat- 
ter halo (Sofue 1998). The extreme options are a potential 
dominated by either a central mass within a radius of ~ 2 kpc 
or a dark matter halo extending out to 100 kpc, the two of 
which crudely limit the escape velocity to lie in the range 
Vesc ^ (1.5 - 3.5) X Vi-ot. For a peak rotational velocity of 
Vi-ot ^ 130 kms"' at a radius of ^ 500 pc (Gotz et al. 1990), 
the escape velocity is Vesc ~ 200 - 460 kms"' , much lower 
than the allowed range of Vqo = 1410 - 2240 km s"' . 

Our estimated terminal velocity for M82's wind is also 
much larger than estimates of the velocity of the Ha emit- 
ting warm ionized medium, which achieves a maximum 
inclination-corrected velocity of vhq ~ 600 - 650 kms"' 
at heights between |z | ~ 350 to 2000 pc from the nucleus 
jMcKeifli et al.1 [1991 Hiopbefl & Bland-Hawthorn 1998). 
This is not surprising, as the C hevalier & Cleggi (il985i) model 
predicts a large difference between the velocities of the wind 
fluid and the velocity of the WIM (which they identified as 
clouds of cooler denser gas entrained and accelerated by the 
wind fluid). In the case of M82 velocity of the WIM com- 
ponent of the wind appears to exceed the galactic escape ve- 
locity, although more generally the WNM and WIM in super- 
winds has a speed more closely comparable to Vesc and hence 
its lon g term fate is ambiguous (iHeckman et alJlioOOHMartinl 
l2005h . 

Our results constrain the energy per particle within the star- 
burst region, neglecting the possibility of radiative energy 
losses and additional mass-loading of the wind fluid exte- 
rior to the starburst region. Is it plausible that these effects 
could be of sufficient magnitude to reduce the terminal ve- 
locity of the wind fluid at large radius to a value below the 
escape velocity? For a starburst with a SN thermalization ef- 
ficiency of e = 0.3 and a central mass loading factor (3=1 the 
larger scale wind would have to lose ~ 90% of its remain- 
ing energy to reduce the final value of Vqo to 460 kms"' in 
the absence of additional mass loading of the wind fluid (98% 
energy loss to reduce v'oo to 200 kms"'). The required en- 
ergy loss fraction is very similar for a central starburst with 
e = 1 and /3 = 2.3. We know that radiati ve energy losses in 
M82's large-scale wind ai-e small, < 10% (iHoopes et al.l2003l: 
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[Strickland et aDl2004bl) . It seems unlikely that the wind fluid 
could so effective in transferring enough of its energy to a 
different gaseous phase, such as the galactic halo medium, 
without a significant fraction of this energy being radiated 
and hence observed, unless the radiation is spread over a re- 
gion much larger than the currently studied 12 kpc extent 
of M82's wind. The remaining alternative is to strongly mass 
load the wind fluid once it has left the starburst region. As- 
suming that energy losses in the wind fluid are negligible then 
in the case of the starburst with e = 0.3 and /3 = 1 we would 
need to effectively mix an additional Mexaa ~ 13 Mq yr"' into 
the wind fluid outside the starburst region to reduce v^o to 
460 km s"' 70 Mq yr"' to reduce Voo to 200 km s"' ). For 
e = 1 and = 2.3 then M^xtra ^ 43 Mq yr"' (voo = 460 km s"' ) 
and 240 Mq yr"' (voo = 200km s"'). Even the smallest of 
these values is uncomfortably large in terms of the required 
gas reservoir, and furthermore the resulting X-ray emissivity 
of the volume filling wind fluid might then exceed the ob- 

. 3 

served soft X-ray luminosity of the wind (EI cx M ). Given 
the preceding arguments we believe that is probable that the 
metal enriched wind fluid will ultimately escape M82. 

5.4. Implications for Understanding Feedback and 
Superwinds 

It is worth stressing that a multi-phase wind model leads 
to phase-dependent ejection efficiencies for the material in 
starburst-driven superwinds. Current cosmological models of 
galaxies that attempt to include superwinds invariably have 
single phase winds (although the simulation may allow mul- 
tiple ISM phases, their galactic winds lack the fine-scale mul- 
tiphase nature of real winds). Cosmological simulations lack 
the ability to separately model the differential ejection of the 
metal-enriched high-energy-per-particle wind fluid from the 
less enriched, lower velocity, but more massive WNM and 
WIM. Their results on metal transport and feedback processes 
in galaxies are potentially misleading or incorrect. 

In particular, simulations that add all the gas mass of the 
WNM and WIM phases to the metal-enriched SN ejecta re- 
quire nonphysically large fractions (e > 100%) of all the avail- 
able SN mechanical energy in star-forming galaxies (normal 
galaxies in addition to starburst galaxies) in order to match 
the observed degree of metal loss from galaxies an d metal 
enrichment of the IG M. (e.^ . Springel & Hernq uis^ l2003bl: 
lOppenheimer & Pavel |2006;KobayashieLaLi |2007V These 
papers interpret such problems as implying the physical inad- 
equacy of purely SN-driven winds and need for addition en- 
ergy or wind sources (e.g. radiation-driven winds). Our inter- 
pretation is that the flaw lies in the inadequacy of their recipe 
for winds. In contrast simulations that effectively treat the 
SN ejecta alone without additional WIMAVNM mass load- 
ing achieve good matches to observed Lya and O VI absorber 
statistics while using reasonabl e SN efficiencies (e ~ 3 0%, 
ICen et al.ll200l ICen & Ostrike3l2C^ ICen & Fang|l2006h . 

We further note that it is unrealistic to apply a single uni- 
form feedback prescription to all star forming galaxies. As 
discussed previously there are good reasons to expect that e 
is a function of the local intensity of star formation. Fur- 
thermore only galaxies with average star formation rates'" 

This value for the critical SsF is based on a Salpeter IMF with mass 
limits of 1 and 100 Mq . 'Dahlem et alj 120061) demonstrate that the star for- 
mation intensity threshold to establish extra-planar gaseous activity depends 
on the host galaxy mass, but as this mass dependence is relatively weak we 
have ignored it in presenting the simple superwind recipe. 



per unit area S sf x'. 4 x 10 ^M^yr ' kpc ^ show evidence 
for superwinds JLehnert & Heckmaniri996bl: lHeckmanll2003l: 
iDahlem et al .l l20()6h. 

In the spirit of iMartinI ( 11999 ^) we now present a simplified 
bimodal feedback and superwind recipe appropriate for use 
in semi-analytical and numerical simulations of galaxy for- 
mation. This recipe is most applicable to simulations that aim 
to answer questions regarding the galaxy mass-metallicity re- 
lationship, and/or the metal content of the IGM, as it primarily 
addresses the metal-enriched stellar wind and SN ejecta from 
starburst-driven superwinds. 

The formulation presented below differs from most cur- 
rently used recipes for winds. It prevents massive galaxies 
with moderate star formation rates (such as the MiUcy-Way) 
and the bulk of the galaxy population at low redshift from 
having any wind-related metal-loss, but the chances of metal 
ejection from intensely star-forming galaxies (starbursts) is 
enhanced. It is unlikely to act as an effective mechanism to 
suppress star formation in all star-forming galaxies, unlike the 
wind recipes in the simulations we criticized earlier 

In galaxies with average star formation rates per units area 
greater than Esf > 4 x lO"^M0yr"' kpc"^ (i.e. starbursts) 
feedback is efficient and leads to the creation of galactic su- 
perwinds. The hot gas in a superwind can be characterized 
by a mass injection rate Mhot of material of metal abundance 
Zhot, with an associated mechanical energy injection rate E^ot- 
In terms of the star formation rate S in units of Mq yr"' for a 
Salpeter IMF between the mass limits of 1 and 100 Mq, then: 

• Mhot ~ 0.3 X /3hot C'5 Mq yr"' , where /3hot is in the range 
1.0-2.8. 

• Zhot is given by Equation [12] and is a function of the 
metal abundance in the merged SN and stellar wind 
eject ZsN+sw, the metal abundance in the cool ambient 
ISM Zcoid and the mass loading factor /3hot- 

• The effective mechanical energy injection rate, ac- 
counting for the efficiency of thermalization, Ehot ~ 
6.5 X 10hQ(C'5ergs"' , where Chot is in the range 0.3 - 
1.0. 

• Alternatively, if winds are implemented by imparting 
a fixed velocity to gas "particles," then the appropriate 
velocity is Vhot ~ 2650(ehot//?hot)'''^ kms"' , and should 
be in the range 1400 - 2200 km s"' . 

• The values of ehot and /3hot are correlated in the sense 
that low values of e are associated with low value s of /3, 
such that Tc or Voo lie in the range described in S I5.1I & 
E3] 

• The participation factor C is the fraction of the total star 
formation that can effectively participate in driving a 
superwind. In the case of M82 C is unlikely to be lower 
than ~ 0.5, and is more plausible in the range ( ^ 0.7 - 
1.0 (see AppendixO. 

In more quiescent star forming galaxies with Ssf < 4 x 
10"^MQyr"' kpc"^ feedback is energetically inefficient and 
no true wind arises, although SN-heated and enriched gas may 
vent into the halo to form a galactic fountain. Galactic foun- 
tains are considerably less constrained theoretically than su- 
perwinds, but it is plausible that all the material eventually 
returns to the galactic disk after some delay. 
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This recipe captures two essential inter-related features of 
starburst and feedback physics that have been established in 
the local Universe and that have not been applied in existing 
recipes for winds in cosmological simulations: the existence 
of a threshold in star formation activity per unit disk area nec- 
essary to drive global galactic winds, and that the quantitative 
properties of feedback are not the same in all galaxies but de- 
pend on the intensity of star formation. 

5.5. Measuring feedback parameters in other starburst 
galaxies 

One limitation of the superwind recipe as presented is 
that it treats all starbursts equally, and applies an M82-like 
wind to them all. Although simulation s of the type made by 
iMelioli & de Gouveia Dal Pinol (|2004 suggest that the effi- 
ciency of supernova feedback may be bimodal, other evidence 
suggests that not all starburst-driven fl ows are equally energy- 
efficient (see e.g. 'Calzet ti et al]l2004. 

The techniques described in [Strickland & H eckmanI (l2007|) 
and in this paper can be applied to other starburst galaxies in 
order to measure their effective supernova thermalization effi- 
ciencies and mass-loading parameters. In general one would 
be limited to a purely spectral X-ray analysis, as for more dis- 
tant starbursts it would be difficult to spatially separate the 
hard diffuse X-ray emission from the emission from X-ray 
binaries and nuclear AGN even with the Chandra X-ray Ob- 
servatory. 

If our model is correct then we expect that the luminosity of 
the diffuse hard X-ray (HX) emission of a starburst region of 
star formation rate S and effective radius 7?^ to scale roughly 
as 

(-2 c2 

Lhxcx£/xAcx^-^— , (16) 

where we have substituted the star formation rate into Equa- 
tion[T5] and have ignored the e and (3 dependence of A. 

This is notable as it predicts Lhx from the wind fluid has a 
stronger than linear dependence on the star formation rate. It 
should therefore be possible to extract this emission compo- 
nent from the other hard X-ray emission components in star 
forming galaxies that have been shown to scale linearly with 
star formation rate (Ranalli et al. 2003; Gri mm et al 1120031 
IColbertetal.]|2004t fPersic & Rephaeli,2007i) . The task wifl 
also be simplified in that it is the hard X-ray line emission 
that has proven itself to provide the strongest constraints on 
e and (3 in M82, rather than the hard X-ray continuum emis- 
sion that the X-ray binaries dominate in the integrated X-ray 
spectra of galaxies. 

We note that Fe Hea emission has been dete cted in the X- 
ray spectra of the nearby starburst NGC 253 ('Weaver et alj 
12002) . the merging system Arp 299 (Balloetal. 2004^ 
and the ultra luminous IR gala xies (ULIRGs) A rp 220 
(llwasawaetal.1 [20051). M rk 273 (Balestraet allllool . NGC 
6240 (lNetzeretal.ll2005h and the Superantennae (Braito et 
al. 2008, submitted). Previously this emission has been inter- 
preted in terms of ion ized reflection associated with an AGN 
jRoss & Fabianll2005 ). and while this is likely true in objects 
with strong AGN (e.g. NGC 6240), ionized reflection might 
not be the cause of the Fe Hea emission in all of these ob- 
jects. These papers often do consider a starburst origin for the 
ionized iron emission, but dismiss it based on the expectation 
that its luminosity should scale linearly with the star forma- 
tion rate. 

Apparently diffuse hard X-ray emission has also been seen 



in the C/iQ«t/ra obser vations of the sta rburst galaxies N GC 
2146 (llnui et al.ll2005h and NGC 3256 (iLira et al.ll20()2h . al- 
though these papers do not discuss hard X-ray line emission. 
The XMM-Ne wton X-ray sp ectrum of NGC 3256 presented 
as Figure 2 of iJenkins et al.l (12004) shows a hint of a hard X- 
ray iron line, but the paper does not quantify the energy and 
intensity of the feature. 

To crudely test the plausibility of a wind fluid origin for the 
Fe Hea emission in these other starbur sts we consider NGC 
253 (LpeHea = (5.3!^|) X lO^crgs'^, IWeaver et al.1 l2002h 
Arp 299 (LpeHea = (9.5 ±4.5) x lO^** erg s"' . iBallo et al.l200l 
Arp 2 20 (LpeHea = (1.2 ±0.6) X 10^° erg s-' Tlwasawa et al.1 
'2005') and Mrk 2 73 (LpeHea = (9.9!^-^) x lO^^^ergs"!, 
Balestra et all l2005h . We wifl use the A = 8 - 1000 fim 
IR luminosities of these galaxies, taken from the data in 
ISanders et al.1 (l2003h . as proxies for the star formation rate 
S, making no correction for the fraction of the IR luminosity 
contributed by any AGN. The LpeHea luminosities given above 
are based on the 6.7 keV line photon fluxes quoted in the ap- 
propriate paper and scaled to the adopted distances given be- 
loW; 

IWeaver et alJ (l2002h found that the hard X-ray emission in 
NGC 253 (D=2.6 Mpc, Lir = 1 .9 x 10'" Lq) comes from a 300 
pc by 60 pc region associated with the nuclear starburst region 
and gas torus. If we interpret this as a disk of that diameter 
and height then the equivalent 1 -dimensional effective radius 
is Ri, = 126 pc. Scaling M82's observed Fe Hea emission by 
(5ngc253/5m82)^ X (/?*.M82/^*,NGC253) we obtain a predicted 
Fe Hea luminosity associated with the starburst of LpeHea ~ 

4.4 X lO-'^ erg s"' , which is consistent with the observed value 
given above. 

Arp 299 (D=48.5 Mpc) is a merging system with a total 
IR luminosity of Lir = 7.7 x 10" Lq, although the Fe Hea 
emission is only associated with the eastern nucleus (IC 694). 
Arp 220 (D=74 Mpc) and Mrk 273 (D=158 Mpc) have com- 
parable IR luminosities of Lir = 1.3 x 10'^ L©. As the sizes 
and geometries of the starburst regions in these galaxies are 
not well established we will ignore the scaling and scale 
the observed M82 Fe Hea luminosity only by the square 
of the ratio of the IR luminosities. For Arp 299 the pre- 
dicted Fe Hea luminosity is LpeHea ^ 3.0 x 10 erg s"' , a few 
times higher than the observed value of ^ lO'^^'ergs"' . For 
Arp 220 and Mrk 273 the predicted luminosity is LpeHea ~ 

8.5 X lO^'^ergs"' . This value is comparable to the (poorly- 
constrained) Fe Hea luminosity of Mrk 273, and several times 
larger than the observed Fe Hea luminosity of Arp 220. As 
we would expect the effective size of the starburst regions in 
these objects to be larger than that of M82 we would expect 
the luminosity to somewhat less than the predictions given 
above, so these predictions are in reasonable agreement with 
the observations. 

Such back of the envelope calculations certainly do not 
prove that the Fe Hea emission seen in these objects arises 
from the same processes as in M82. They do demonstrate 
that there is some hope of being able to test and explore this 
model in objects other than M82 with current X-ray data from 
the Chandra and XMM-Newton X-ray Observatories. 

In the longer term the superior spectral resolution and sensi- 
tivity of the calorimeter instrument on the International X-ray 
Observatory will allow routine detection of hard X-ray line 
emission from a much larger sample of star forming galax- 
ies. Furthermore, with an energy resolution of AL 4 eV at 
E = 6.7 keV it will allow us to test the assumption that these 
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lines arise in coUisionally heated gas. 

6. CONCLUSIONS 

We have measured the net energy efficiency of supernova 
and stellar wind feedback (the thermalization efficiency e) 
in the starburst galaxy M82 and the degree of central mass- 
loading j3 of its superwind by comparing a large suite of 1 
and 2-dimensional hydrodynamical models to a set of obser- 
vational constraints derived from hard X-ray observations of 
the starburst region (the fluxes of the Hea and Lya-like lines 
of S, Ar, Ca and Fe, along with the total diffuse £ = 2 - 8 keV 
X-ray luminosity). 

The observational constraints are based on a combination 
of the best currently available observations of the central 
500 pc of M82 with both the Chandra and XMM-Newton 
X-ray Observatories. We com bine the results presented in 
[Strickland & HeckmanI (l2007l) on the diffuse hai'd X -ray Fe 
line emission and the diffuse £ = 2 - 8 keV hard X-ray lumi- 
nosity with new results (this paper) on the S, Ar and Ca lines. 
By simultaneously applying a large number of observational 
constraints (we used eight different line fluxes or upper lim- 
its, an upper limit on the broad-band X-ray luminosity and a 
physical requirement that the flow not be strongly radiative) 
to the simulated starbursts we have been able to constrain the 
physical properties of the starburst region more tightly than 
by using more standard X-ray diagnostics. 

We modeled the physical properties of the SN-heated and 
enriched plasma within an assumed disk-like starburst re- 
gion and the flow of gas outside this region out to a ra- 
dius of 500 pc usin g either an appropr iately scaled version 
of the 1-dimensional lChevalier & Cleggl(il98 5) analytical su- 
perwind solution, or with a full hydrodynamical simulation 
using the VH-1 hydrodynamics code in 2-dimensional cylin- 
drical symmetry. We demonstrated a method for correcting 
the 1 -dimensional radial CC model so as to apply to a disk- 
like starburst geometry with reasonable accuracy. Although 
this introduces errors of < 30% in the individual predicted 
hard X-ray line fluxes compared to the full 2-D hydrodynam- 
ical solution, such errors have little effect on the derived con- 
straints on e and (3 as the uncertainties in the observed line 
fluxes can be larger than this. We also considered the issue 
of the applicability of the adiabatic CC wind solution and find 
that for the study of the volume-filling gaseous component of 
the M82 superwind within individual star clusters and over 
the starburst region as a whole the use of adiabatic models is 
valid. 

We considered a broad range of plausible parameters for 
the M82 starburst, varying the age and mode of star forma- 
tion, the starburst region size and geometry, and supernova 
metal yields. We further considered a correction accounting 
for the fraction of massive stars that occur in high density re- 
gions where high radiative energy losses will preclude their 
making any contribution to the superwind, which we term the 
participation fraction Q- 

We find that over all these varied input parameters (we 
calculated 663 full hydrodynamical simulations using differ- 
ent model parameters, and 4615 models using the scaled CC 
model), those models that do satisfy the existing observational 
constraints all have medium to high SNh-SW thermalization 
efficiencies (30% < e < 100%) and have mildly centrally 
mass-loaded wind fluids (generally 1.5 < /3 < 2.5, although 
in a few cases /? = 1 is allowed). 



Our results are ultimately based on the measured Fe 
Hea emission of the SN heated and enriched plasma driving 
M82's wind. As such they represent the first direct measure- 
ment of the efficiency of SN feedback and mass-loading of 
the wind fluid in a starburst galaxy. Future observations with 
higher spectral resolution (e.g. with ASTRO-H and IXO) will 
be required to test whether the observed line emission truly 
arises from a collisional plasma in or close to ionization equi- 
librium. 

That the majority of the models that successfully match 
our imposed hard X-ray observational constraints also predict 
thermal pressures in the starburst region (7.0 < logP^/k < 
7.5) that are consistent with prior independent estimates of 
the starburst region pressure increases our confidence in the 
validity of our method and results. 

The allowed values of the SN and SW thermalization effi- 
ciency and the mass loading factor correspond to a tempera- 
ture for the metal-enriched plasma within the starburst region 
of 30 - 80 million Kelvin, a mass flow rate of the wind fluid 
out of the starburst region of Mtot ^ 1.4 - 3.6 M0 yr~' and 
imply a terminal velocity of the wind (in free expansion) in 
the range Vqo = 1410 - 2240km s"' . This velocity is con- 
siderably larger than the gravitational escape velocity from 
M82 (vesc < 460 km s"') and the velocity of the Ha emit- 
ting clumps and filaments in the inner 2 kpc of M82's wind 
(vHa ^ 600 km s"' ). Reducing the velocity of the wind fluid 
once outside the starburst region to the escape velocity would 
require either unrealistically large mass loading (Mextra '^-^ 13 
- 240 Mq yr"') or excessive energy losses (^ 90 - 98%) that 
must also avoid raising the luminosity of the wind. Thus the 
hot gas driving M82's wind will probably escape the galaxy's 
potential well. This result also adds support to the idea that the 
long term fate of gas within superwinds is phase-dependent. 

We show that the luminosity of the hard X-ray line emission 
produced by the wind fluid in a starburst region scales with the 
star formation rate S and starburst size R^, as Lux S^/R^,. 
This should serve to distinguish this emission from the other 
sources of both soft and hard X-ray emission in star forming 
galaxies that scale linearly with the star formation rate. We 
show that existing measurements of Fe Hea emission from 
other actively actively star-forming galaxies (NGC 253, Arp 
299, Arp 220, Mrk 273), previously interpreted as ionized re- 
flection from their AGN, are also consistent with our predic- 
tion for the wind fluid given the existing observational uncer- 
tainties. 

Drawing on our results for M82 we have also presented a 
simple prescription for implementing starburst-driven super- 
winds in cosmological models of galaxy formation and evolu- 
tion that more accurately represents the energetics of the hot 
metal-enriched phases than do existing recipes in the litera- 
ture. 
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APPENDIX 

A. THE CONDITIONS UNDER WHICH THE ADIABATIC SOLUTION BECOMES INVALID 

The basic adiabatic CC model is a function of only three parameters: £tot, Mtou and R^,. ISilich et al] (l2004t) were the first to 
point out that for certain values of these parameters, in particular small 7?* and/or high Mtot, the total luminosity predicted by 
the CC model was comparable to the net energy injection rate and that the analytical CC solution was no longer self-consistent 
as it assumed energy losses to be negligible. They de veloped a self-consistent radiative solution, later tested with numerical 
hydrodynamical simulations (" Tenorio-Tagle et ani2007h . 

There in not an exact and easily evaluated function that may be used to assess whether the adiabatic assumption is valid for any 
particular combination of Etou Mtou and Nevertheless. ISilich et alJ (l2004l) state that the threshold between the adiabatic and 
the fully radiati ve solution occurs when the radiative energy losses per unit volume exceed ^ 30% of the energy injection rate 
per unit volume. iTenorio-Tagle et all (l2007l) show that this is most likely to occur at the very center of the star cluster or starburst 
region. 

In terms of the variables defined in § |3] the adiabatic solution becomes invalid when 

,ZwF /\^'-^^tOt /A1\ 

2^4 ^ y* 



Equivalently, if 



then radiative losses are high and the CC model is invalid. Note the strong dependence on the mass injection rate and any mass 
loading, and inverse dependence on the energy injection rate, thermalization efficiency and effective radius. The higher the mean 
energy per particle, or the larger the starbursting region is, the less-likely significant radiative losses become. 

A. 1 . Individual star clusters 

ISiUch et al] (|2004|) and lTenorio-Tagle et al.l (l2007h provide examples of star clusters with properties such that radiative losses 
are severe. However, the mass and energy injection rates they generally use have significantly lower values of the energy per 
particle (and hence higher radiative losses for a given assumed metal abundance) than what we consider appropriate for clusters 
in a M82-like starburst. 

Equation IA2I mav be rewritten as an explicit function of the cluster mass. We rewrite the energy and mass injection rates 
as a function of the star cluster mass assuming in this case an instantaneous starburst model, such that Esn+sw = jM(, and 
^SN+sw = kM(y, where M(y is the total cluster mass in units of 10^ Mq (again using a Salpeter IMF between the lower and upper 
mass limits of 1 - 100 M©). Note that the value of j and k are a functions of the age of the cluster, and are largest in magnitude for 
a cluster age of f ~ 4 Myr in the Starburst99 models. The critical cluster mass above which cooling exceeds 30% of the energy 
injection and the CC model is no longer valid is then 

^h4/^.c1 

0.74Cd/5e'ifc3A(r,ZwF) 

For the purposes of illustration we consider a cluster of age f = 5 Myr, for which jM(, = 3.7 x 10'*" ergs"' and kM(, = 2.3 x 
10"^ M0 yr"'. For unit e, /?, and C the central temperature in the CC model would be ^ 7.5 x lO^K (equivalent to v'oo = 
2250km s"' ). Assuming Zsn+sw = 5Zq the cooling function is A(r, ZwF) = 5.54 x lO"^-' erg cm-' s"' . 

The massive star clusters observed by iMcCrady & GrahamI (l2007h in the nuclear starburst region of M82 have half-light radii 
of between 1 - 3 pc. We adopt a typical cluster radius of iJ^ci = 2 pc. Substituting these values into Equation IA3 1 we find 

,2 4i / 5.54xl0-" \ fR^ 

CciPi V A(r,ZwF) y V2pc 

For unit thermalization efficiency and no additional mass loading then no realistic cluster will lie in the radiative regime. Signifi- 
cant mass loading within the cluster (e.g. f3ci = 5) or low thermalization efficiency (ecI = 0.1) are required to reduce M(y to of order 
unity for a f = 5 Myr old cluster. 

An older cluster is even less likely to be radiative. From the Starburst99 models we find that at an age of 20 Myr the mass 
injection from stellar winds has declined dramatically yielding a net mass injection rate a factor 5 lower than at an age of 5 Myr 
The energy injection rate is also lower, but only by a 2.5. At this age even clusters of mass a few times 10^ Mq are still adiabatic 
even if /3ci = 5 or ed = 0. 1 . 

A. 2. Starburst regions as a whole 

We rewrite the energy and mass injection rates as a function of the star formation rate S, so that /isN+sw = oS and Msn+sw = bS. 
If we use the values given in Table H] for constant star formation with an age of 30 Myr, again assuming a metal abundance of 
ZsN+sw = 5 Z0, we find that the star formation rate must exceed 



Me > 1.9 X 10^-4^ r.; " ■ (A4) 
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Fig. 12. — Gas number density and velocity vectors in two hydrodynamical simulations of a starburst in (a) a spherical starburst of radius = 300 pc and (b) a 
disk-like starburst with radius r» = 375 pc and total thickness ht = 105 pc. The central gas density and pressure is the same for these configurations. The highest 
density contour shown corresponds to logn = —1.00 and further contours correspond to successive decreases in log« in steps of 0.25. 

before the adiabatic CC solution becomes invalid. 

For M82 the star formation rate is 5 4 - 6 Mq yr"' for the Salpeter IMF assumed, and the starburst region has effective radius 
of ^ 100 pc at smallest (and more likely ~ 300 pc). Thus the CC model is valid unless mass loading is significant /? > 5) 
or thermalization is inefficient e < 0.1, or some combination of the two. Note that a low value for the participation fraction C, will 
reduce the tendency of the starburst region to be radiative. 

We conclude that under most, but not all, conditions the adiabatic IChevaher & Cleg3 (119851) model is appropriate for M82's 
nuclear starburst. 

B. CORRECTIONS TO THE CC MODEL DUE TO A NON-SPHERICAL STARBURST REGION 

It is possible to use the spherical analytical CC model to predict the wind-fluid properties and emission from a disk-like starburst 
of diameter c/* and total height /i* by applying modest numerical corrections that are a fixed function of the geometry alone. As 
long as the adiabatic assumption is valid these correction factors do not depend on the starburst energy or mass injection rates, 
and are thus independent of the thermalization efficiency e and the mass loading factor /3. 

In the analytical CC model the surface area of the starburst region S^, controls the gas density with in the starbu rst region by 
limiting the rate at which matter can flow out of it, as described in § 13.21 (see in particular Equ. 30 of ICanto et al.ll20 00). Note 
that other parameters, in particular the central temperature and the terminal velocity of the wind are independent of the starburst 
geometry. Thus, for a given fixed energy and mass injection rate, the central density and pressure in a disk-like starburst of 
surface area 5* = 7r(i*((i*/2 + /z*) should be very similar or identical to values in a spherical starburst of the same surface area 5*. 

We tested whether this is indeed the case using the VH-1 hydrodynamical code run in 2-D cylindrical symmetry as described in 
§ 13.11 This code reproduces the standard spherical CC solution to within an accuracy of 0.1% for an assumed spherical starburst. 
For a disk-like starburst of diameter d^, = 750 pc and = 105 pc (§ |3.3.2| l a spherical starburst of equivalent surface area has 
= 300 pc. 

We find that the central density and pressure predicted by VH-1 for this disk-like starburst geometry differs from the values 
it predicts for a = 300 pc spherical starburst by only 3.1%. Adopting another disk-like geometry with the same surface area 
(c/* = 682 pc and /i* = 200 pc) yields central density and pressures within 4.2% of the spherical solution. For the other disk-like 
starburst geometry suggested by observations of M82, with di, = 1015 pc and hi, = 140 pc the 2-D simulations predict central 
densities and pressures within 3.7% of the spherical starburst with equivalent surface area (/?* = 405 pc). These simulations 
demonstrate that we can accurately predict the gas properties within a disk-like starburst region using the spherical CC model of 
equivalent starburst surface area. 

The flow pattern exterior to the starburst region is also influenced by the geometry of the starburst region. A disk with c/* > 
will weakly collimate the flow along the minor axis and give the superwind a preferred direction even in the absence of external 
coUimation by the ISM. In Fig. [12] we illustrate this by plotting the density and velocity fields for equivalent spherical and 
disk-like starbursts. 

For our purposes the main effect is that the resulting non-spherical flow alters the total volume emission integral within any 
fixed region that is larger than the starburst region. Within the starburst region itself this effect is not particularly significant when 
evaluating total gas masses or emission integral given the relative uniformity of the density and pressure. However, outside and 
above a disk-like starburst region the gas density drops off with increasing height more slowly that the p cx R~^ predicted by the 
CC model. 

Calculating the correct emission integral for a disk-like starburst based the results of the spherical CC model thus entails some 
geometry-based corrections because of the different starburst region volumes and secondly due to the non-spherical flow pattern 
outside the starburst region. The first correction is trivial to evaluate analytically. The second correction must be determined from 
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a full hydrodynamical simulation of the specific disk-like starburst, and further depends on the size of region Robs we wish to 
evaluate the emission integral over. 

Our procedure for calculating the gas properties and total emission of disk-like starburst regions using the spherical CC model 
is as follows: 

1 . For a given starburst geometry evaluate the surface area S^„ and run a simulation using the CC model using the effective 
starburst radius = (5'*/47r)'/^. This accurately predicts the values of fluid variables such as density, temperature and 
pressure within the starburst region. 

2. Macroscopic variables, in particular the volume emission integral, require further geometrical corrections to make the emis- 
sion integral predicted using the spherical CC model, EIqc, match the emission integral found in a full multi-dimensional 
simulation of a non-spherical starburst region fi/nonsph- 

3. The volume of the starburst region in the spherical model is larger than the volume of the disk-like starburst region, so the 
volume emission integral predicted by the CC model within some radius Robs > Ri, must be scaled by a factor Ky equal to 
the ratio of the starburst region volumes: 

Ky= (Bl) 

4. A full hydrodynamical simulation for a given disk-like starburst must be run to evaluate the correction /Tnonsph = 
EInons^\i/{KvEIcc) to the emission integral within the chosen Robs due to the non-spherical flow pattern manifest out- 
side the starburst region. We use Robs = 500 pc to match the region from which the M82 nuclear region spectra analyzed in 
Paper I were extracted. Simulations of different disk-like starburst regions where we altered the thermalization efficiency 
and mass loading demonstrated that the correction /Tnonsph to the total emission integral is a fixed function of the geometry, 
and is independent of the energy and mass injection rates. 

5. The temperature-dependent differential emission measure for the scaled CC model EIscciT) is calculated using the factors 
discussed above: EIscciT) = Ky x A'nonsph x Ef^'-iJ). Combining this with the temperature-dependent cooling functions 
and summing over all temperatures gives the final predicted luminosity for the chosen energy band or line. 

For Robs = 500 pc and a disk-like starburst with d^, = 750 pc and = 105 pc, the primary geometry we apply to M82, the 
non-spherical emission integral correction factor Anonsph = 1-64. For the alternate disk-like starburst geometry with dj, = 1015 pc 
and h-i, = 140 pc the correction factor is A'nonsph = 1.52. For the compact starbursting disk with = 300 pc and hi, = 60 pc the 
correction factor is /Tnonsph = 1 -59. 

With these scalings the spherical CC model can be used to predict to high accuracy the total emission integral from the 
wind fluid within Robs for a disk-like starburst. This method is less accurate at reproducing the luminosity predicted by a full 
2-dimensional hydrodynamical simulation because non-spherical flow pattern leads to subtle differences in the shape of the 
differential emission integral EI(T) compared to a spherical CC solution. 

We compared the luminosities predicted by a small suite of hydrodynamical models with different disk- like starburst geometries 
to the results from the scaled CC analytical model in order to assess the significance of this effect. The error introduced by the 
effect is small in terms of computing the total cooling rate (< 2% for the Sutherland & Dopita cooling curve), or even the the 
X-ray luminosity over the broad 0.3 - 2.0 keV or 2 - 8 keV energy ranges (< 10%). For the high energy lines we are interested in 
because of their temperature sensitivity the effect is larger, and the associated errors in predicted line luminosity can be as large 
as - 30%. 

Nevertheless these uncertainties do not compromise the usefulness of the scaled CC model as they are comparable to the 
observational uncertainties in the nuclear region S, Ar and Ca luminosities of the well detected lines from the XMM-Newton data, 
and less than the uncertainties in the diffuse Fe Hea line luminosity derived from the Chandra ACIS-S data. The relative effect 
of these uncertainties on the region of model parameter space matching our chosen constraints can be seen in Fig.|2] 



C. ESTIMATING THE PARTICIPATION FRACTION C 

There are two main environments that in principle might reduce the fraction C, of all the SNe and stellar winds in a starburst 
region that can contribute to powering a superwind. The first, likely to to be the most significant but also hardest to estimate, 
is physical confinement and increased radiative cooling of young supernova remnants and stellar winds occurring in regions of 
high ambient gas density such as molecular clouds. The second is that radiative losses may be high for those SNe and stellar 
winds associated with massive stars in massive and compact star clusters (particularly super star clusters), where the adiabatic 
CC solution is no longer valid. 

There are some observational findings that ar e suggest that some fraction of the SNe in M82's nuclear starburst occurred withi n 
regions of gas number density n > 10^ cm"^ jKronbe rg et al. 2000; Chevali er & Franssonll200l1; ISeaquist & Stankovicll200'7l) . 
although there are problems with this interpretation (see .Beswick et al..,200& and references therein). However, based on the 
total molecular gas mass of ^ 10^ Mq wit hin the starburst region, the fraction o f the total starburst region volume occupied by 
such dense gas must be 77 < 0.1 or lower (' Lugten et al.lll986l : iNakai et alj|1987l) . Massive stars will be born in regions of high 
gas density, so that if confinement does occur the magnitude of the drop in ^ may well be larger the the filling factor of the dense 
gas. At present we can not rigorously quantify the effectiveness of confinement, nor the fraction of stars that might be affected. 
IChevalier & FranssonI (1200 1) speculate that perhaps as much as 50% of the SNe in M82 occur in dense regions, in which case 
the participation factor may be as low as C 0.5. 
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It is possible to quantitatively estimate the role that over-cooling of the thermalized SN and stellar wind ejecta in dense star 
clusters have on (. Based on the discussion in Appendix lAl we conclude that only the compact clusters of radii a few parsecs 
and masses Miad 10^ Mq (Salpeter IMF, mass limits 1 - 100 Mq) are likely to suffer significant radiative losses (and then only 
under certain conditions). 

For a star clu ster mass function of the form dNc\ oc dM between the mass limits Md.min and Ma.max and where 7 = 2 
(IWeidner & Kroupa.2006,) . the fraction of the total stellar mass in clusters with Md > Miad is 

InMcl.max-lnMrad ^^^^ 
lnMd,max-lnMmin' 

We adopt a minimum cluster mass of Ma.mm = 5 Mq and a maximum cluster mass of Mei,max = 10^ dWeidner et al.ll2004 . 
For Mrad = 10^ Mq then / = 0. 16. Thus the magnitude of the reduction in ( due to radiative cluster winds is small, of order 15%, 
even if SN and stellar wind thermalization is inefficient, and/or there are long-lived sources of additional mass loading in the 
massive star clusters. 
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